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Abstract 

This  study  investigated  the  use  of  ballistic  and  “Double-Dip”  aerobraking  reen¬ 
try  to  return  the  Space  Maneuver  Vehicle  (SMV)  from  geotransfer  orbit  in  no  more 
than  two  atmosphere  passes.  Lift  and  drag  accelerations  were  applied  to  the  two-body 
problem  when  either  of  their  magnitudes  exceeded  g.  Lift  and  drag  coefficients, 
along  with  the  SMV  model,  were  taken  from  Investigation  of  Atmospheric  Reentry  for 
the  Space  Maneuver  Vehicle  by  Captain  McNabb,  AFIT/GA/ENY/04-M03.  Target 
perigees  were  formulated  using  the  two-body  problem.  The  orbit  from  each  target 
perigee  was  numerically  integrated  around  a  planar  earth  model  using  a  fourth  order 
Runge-Kutta  method.  Ballistic  and  “Double-Dip”  reentry  schemes  were  attempted 
with  45  and  70  km  altitude  floors.  Ballistic  reentry  produced  a  near  circular,  low 
earth  orbit  when  the  SMV’s  perigee  altitude  resided  between  66.801  and  68.449  km 
for  a  one  pass  reentry  and  72.226  and  73.445  km  for  a  two  pass  reentry.  “Double-Dip” 
reentry  produced  a  near  circular,  low  earth  orbit  when  the  SMV’s  perigee  altitude 
rested  between  62.416  and  64.962  km.  The  resulting  perigee  windows,  their  respective 
heating  rates,  and  experienced  accelerations  were  analyzed.  Effects  of  uncertainty  in 
the  atmosphere  model  on  successful  perigee  windows  for  each  reentry  scheme  were 
analyzed  by  repeating  the  simulation  with  an  increased  atmospheric  density. 
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Investigation  Of  Aerobraking  To  Return  The  Space 
Maneuver  Vehicle  To  Low  Earth  Orbit  From 
Geotransfer  Orbit 

I.  Introduction 


Background 

The  1999  Department  of  Defense  Space  Policy  states, 

Space  is  a  medium  like  the  land,  sea,  and  air  within  which  military  activ¬ 
ities  will  be  conducted  to  achieve  U.S.  national  security  objectives.  The 
ability  to  access  and  utilize  space  is  a  vital  national  interest  because  many 
of  the  activities  conducted  in  the  medium  are  critical  to  U.S.  national  se¬ 
curity  and  economic  well-being.  The  globally  interdependent  information- 
and  knowledge-based  economy  as  well  as  information-based  military  oper¬ 
ations  make  the  information  lines  of  communication  to,  in,  through,  and 
from  space  essential  to  the  exercise  of  U.S.  power.  [3,  2] 

The  United  States’  need  to  effectively  control  space  and  effectively  defend  its  space 
assets  increases  every  day.  Not  only  have  our  war  fighters  become  exceedingly  depen¬ 
dent  on  satellites  for  top-of-the-line  battlefield  awareness  and  communication,  U.S. 
fighters  and  bombers  also  rely  on  the  Global  Positioning  System  (GPS)  constellation 
to  guide  their  smart  bombs,  and  many  of  the  U.S.  intelligence  services  rely  on  satel¬ 
lite  imagery  to  investigate  the  actions  of  countries  deep  within  their  borders.  When 
a  satellite  breaks,  launching  a  new  payload  into  space  can  take  years  following  the 
current  space  launch  paradigm,  but  if  the  U.S.  is  to  effectively  defend  its  space  su¬ 
periority,  a  quick  launch  vehicle  must  be  developed.  Air  Force  Space  Command  has 
suggested  developing  the  Military  Space  Plane  (MSP)  in  order  to  fill  this  need.  The 
MSP  has  four  primary  components  [13,  1-2]: 

•  Space  Operations  Vehicle  (SOV):  a  reusable  launch  element  to  reduce  the  cost 
of,  improve  the  flexibility  of,  and  increase  the  responsiveness  of  earth-to-orbit 
operations; 
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•  Space  Maneuver  Vehicle  (SMV):  an  uncrewed,  reusable  satellite  bus  and  upper 
stage  with  significant  maneuvering  capability; 

•  Modular  Insertion  Stage  (MIS):  a  low  cost  expendable  upper  stage  to  reduce 
the  cost  and  improve  the  operability  of  Space  Support  missions; 

•  Common  Aero  Vehicle  (CAV):  an  aeroshcll  designed  to  deliver  material  (e.g. 
munitions,  UAVs,  or  critical  supplies)  through  space  directly  to  a  theater. 

“The  military  needs  reliable,  responsive,  cost-effective,  and  assured  access  to 
space”  [14,  1],  The  SMV  is  of  primary  concern  in  this  study  because  it  fulfills  the 
military’s  need  stated  above.  “The  SMV  system  concept  combines  the  best  of  aircraft 
and  space  elements  into  a  low  cost,  low  risk  package  which  specifically  addresses 
multiple  warfighter  needs.  With  this  capability,  the  SMV  can  provide  the  United 
States  with  a  flexible  space  superiority  capability  sufficient  to  maintain  battlefield 
dominance  in  support  of  future  warfighter  operations”  [10,  2],  An  operational  SMV 
might  include  [1]: 

•  Up  to  1,200  pounds  of  sensors/payload 

•  72-hours  or  less  turnaround  time  between  missions 

•  Up  to  12-month,  on-orbit  mission  duration 

•  Rapid  recall  from  orbit 

•  Up  to  3.2  km/s  on-orbit  velocity  change  for  maneuvering 

Protecting  and  repairing  friendly  assets,  identifying  and  destroying  space  threats, 
and  placing  new  satellites  into  orbit  are  the  five  most  important  pillars  to  truly  com¬ 
manding  control  of  space.  According  to  AFRL  the  SMV  could  perform  missions  such 
as  [1]: 


•  Intelligence,  surveillance  and  reconnaissance  of  ground  targets  (with  either  in¬ 
tegrated  or  deployable  ISR  satellite) 

•  Deployment  and  recovery  of  microsatellites  (e.g.  Space  Control  Satellites) 
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•  Rapid  constellation  replenishment 

Past  Research 

There  have  been  several  spacecraft  and  concepts  developed  that  have  furthered 
the  technology  needed  to  field  the  SMV.  In  2001,  the  National  Aeronautics  and  Space 
Administration  (NASA)  withdrew  funding  from  the  X-33,  a  single  stage  to  orbit 
vehicle,  and  the  X-34,  a  reusable,  unmanned  suborbital  spacecraft  [4].  The  X-34  flight 
test  model,  built  by  Orbital  Sciences  Corporation,  was  58.3  feet  long,  had  a  wingspan 
of  57.7  feet,  and  stood  11.5  feet  tall.  It  was  to  be  powered  by  a  reusable  Fastrac 
engine,  an  engine  that  provides  60,000  pounds  of  thrust  and  is  aimed  at  boosting 
payloads  up  to  500  pounds  [17].  The  X-34  was  designed  “to  be  air-launched  from 
Orbital  Science’s  L-1011  aircraft,  then  accelerated  to  speeds  up  to  Mach  8,  reaching 
altitudes  up  to  250,000  feet.  It  was  to  land  horizontally  on  a  conventional  runway”  [8]. 
Figure  1  shows  the  X-34. 


Figure  1:  X-34  [17] 


The  X-33  was  envisioned  to  be  69  feet  long  and  have  a  wingspan  of  77  feet. 
Two  J-2S  Linear  Aerospike  engines  were  to  provide  410,000  pounds  of  thrust  to  lift, 
the  285,000-pound  vehicle  to  low  earth  orbit.  The  X-33  was  to  launch  vertically  like  a 
rocket,  reaching  a  suborbital  altitude  at  speeds  greater  than  Mach  13,  and  then  land 
like  an  airplane  [18].  The  X-33  was  supposed  to  be  a  testbed  for  a  reusable  cryogenic 
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tank  system,  cryogenic  insulation,  and  an  integrated  thermal  protection  system  [7]. 
It  was  also  going  to  verify  the  thermal  protection  system  durability,  demonstrate 
guidance,  navigation,  and  control  systems  including  autonomous  flight  control  during 
takeoff,  ascent,  flight,  reentry,  and  landing  [7].  Figure  2  shows  an  artist’s  conception 
of  the  X-33. 


Figure  2:  Artist’s  Conception  of  X-33  [18] 


Although  the  X-33  and  X-34  programs  ended  in  2001,  two  new  programs,  the  X- 
40A  and  X-37,  made  and  are  making  gigantic  leaps  forward  in  automated  spacecraft 
technology.  The  X-40A  is  an  eighty-five  percent  scale  model  of  the  X-37  that  com¬ 
pleted  its  seventh  and  final  automated  landing  test  on  17  May  2001  [6].  The  X-40A, 
shown  in  Figure  3,  successfully  validated  the  low  speed  atmospheric  flight  dynam¬ 
ics  of  the  X-37  at  NASA  Dryden  Research  Center  after  it  was  dropped  seven  times 
from  an  Army  Chinook  helicopter,  acquired  the  runway,  and  performed  a  successful 
landing  [6]. 

The  X-37,  shown  in  Figure  4,  is  being  assembled  by  Boeing  at  their  High  Desert 
Assembly,  Integration  and  Test  Facility  in  Palmdale,  California  [5].  “The  re-usable 
space  plane  will  demonstrate  41  airframe,  propulsion  and  operations  technologies 
aimed  at  significantly  cutting  the  cost  of  space  flight.  The  X-37  can  be  carried  into 
orbit  by  the  Space  Shuttle  or  be  launched  by  an  expendable  rocket”  [21], 

Richard  Powell  led  a  study  for  NASA  in  1985  concerning  an  aerobraking  pro¬ 
cedure  to  return  a  spacecraft  from  geosynchronous  orbit  to  low  earth  orbit  in  order 
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X-40A  Free  Flghl  95 


Figure  3:  X-40A  [6] 


NASA  Dry  den  Flight  Research  ('enter  Photo  Col  lection 
http://www.dlrc.naM.gov/trallcn7Photo/inilcx.litml 
NASA  Photo:  EC99-45I4S-2  Date:  IW 


Illustration  of  the  X-37  Advanced  Technology  Denionstnilor  during  flight 


Figure  4:  X-37  [16] 

to  rendezvous  with  the  Space  Shuttle.  The  Aero-assisted  Orbital  Transfer  Vehicle 
(AOTV)  was  designed  to  leave  geosynchronous  orbit  by  performing  a  small  burn  to 
lower  its  perigee  into  the  atmosphere.  After  aerobraking  during  its  first  dip  into  the 
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atmosphere,  the  AOTV  would  enter  a  300  nautical  mile  phasing  orbit  where  it  would 


perform  a  second  burn  to  lower  itself  into  the  Space  Shuttle’s  orbit.  It  would  circular¬ 


ize  its  orbit  with  a  third  burn  and  rendezvous  with  the  Space  Shuttle  [15,  46].  Figure  5 


shows  the  conceptual  high  lift-to-drag  AOTV,  and  Figure  6  shows  the  AOTV’s  envi¬ 


sioned  mission  profile. 


AOTV  CONFIGURATION  FOR  HIGH  LIFT/DRAG 


Figure  5:  AOTV  Configuration  for  High  Lift/Drag  [15,  44] 


GEOSYNCHRONOUS  ORBIT  MISSION  USING  AN 
AEROASSISTED  TRANS!  iR  VEHICLE 


Figure  6:  Potential  AOTV  Aerobreaking  Flight  Path  [15,  46] 


The  NASA  study  led  by  Powell  lends  a  great  deal  of  plausibility  to  this  study’s 
search  for  an  aerobraking  orbit  that  returns  the  SMV  to  a  low  earth  orbit.  The 
aerodynamic  performance  charts  reveal  that  the  SMV  and  AOTV  share  similar  char¬ 
acteristics  [15,  45].  At  a  perigee  of  70  kilometers,  the  SMV  has  a  coefficient  of  lift  of 
approximately  0.51,  a  lift-to-drag  ratio  around  .807,  and  travels  at  roughly  93  percent 
of  escape  speed.  Comparing  the  SMV’s  traits  to  the  AOTV’s  traits  reveals  that  the 
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AOTV  and  SMV  are  not  identical  but  are  very  similar.  Figure  7  shows  the  AOTV’s 
aerodynamic  properties. 


HIGH  L/D  PERFORMANCE  AERODYNAMICS 


0  10  20  30  HO 

a,  degrees 

Figure  7:  High  L/D  AOTV  Performance  Aerodynamics  [15,  45] 

The  AOTV’s  altitude  history  will  be  used  to  help  determine  what  altitude  to 
start  the  simulation.  While  the  AOTV  was  aerobraking  it  dipped  to  an  altitude  near 
54.8  kilometers  [15,  47]  as  seen  in  Figure  8.  This  study  will  begin  at  45  kilometers  and 
gradually  increase  the  SMV’s  perigee  altitude.  There  was  a  difference  in  Powell’s  study 
that  might  affect  the  outcome  of  this  two  dimensional  study.  Powell’s  simulation  was 
three  dimensional  and  accounted  for  any  inclination  change  required  to  rendezvous 
with  the  shuttle.  The  extra  energy  bled  off  during  an  inclination  change  in  the  Powell 
study,  but  not  accounted  for  in  this  study,  might  lead  to  failure. 

A  large  amount  of  information  is  drawn  from  Captain  Dennis  McNabb’s  thesis, 
ENY/04-M03,  entitled  Investigation  of  Atmospheric  Reentry  for  the  Space  Maneuver 
Vehicle.  His  results  will  be  discussed  in  Chapter  II. 

Concept  of  Operations 

The  current  space  launch  paradigm  prohibits  swift  launch  or  retrieval  missions. 
The  inability  to  quickly  repair  or  replace  friendly  spacecraft  and  assess  or  destroy 
unfriendly  spacecraft  severely  limits  the  Department  of  Defense’s  ability  to  protect 
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ALTITUDE  HISTORIES  FOR  MAXIMUM  RETURN  WEIGHT  AOTV’S 
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Figure  8:  Altitude  History  for  Maximum  Return  Weight  AOTV  [15,  47] 

the  United  States.  Without  guaranteed  space  access,  troops  in  the  field  might  lose 
realtime  communications,  and  hostile  activities  might  go  undetected.  Analysis  of  the 
SMV’s  proposed  capabilities  will  reveal  that  the  SMV  would  fill  that  gap. 

The  SMV  will  be  delivered  to  a  circular,  low  earth  orbit  by  either  an  Evolved 
Expendable  Launch  Vehicle  (EELV)  or  one  of  the  space  shuttles.  Equation  1  gives 
the  required  AV  for  a  coplanar  burn  from  a  circular,  low  earth  orbit  to  geotransfer 
orbit  [22,  305-306], 

AV'  =  ]/2r+S  ~  Vt  (1) 

where 

H  =  398600.5  ( km3/s 2) 
r  =  Magnitude  of  SMV  position  vector  (km) 

£  =  Energy  of  transfer  orbit 

From  a  400  km  circular  orbit,  a  AV  of  approximately  2.4  km/s  is  required  to  enter 
geotransfer  orbit.  Assuming  an  instantaneous  burn,  the  Hohmann  transfer  orbit  would 
deliver  the  SMV  to  the  target  position  or  satellite  in  slightly  under  5.3  hours.  After 
it  arrives  at  its  target  the  SMV  could  perform  any  number  of  tasks.  Once  its  task 
was  accomplished,  the  SMV  would  be  faced  with  the  difficult  task  of  returning  to 
earth  for  any  necessary  changes  or  repairs.  The  SMV’s  advertised  AV  capability  of 
3.2  km/s  is  1.6  km/s  short  of  the  4.8  km/s  required  for  a  two  burn  round  trip.  If 


the  SMV  is  going  to  operate  in  geotransfer  orbits,  an  alternate  return  method  must 
be  developed.  Aerobraking,  the  process  of  dipping  a  spacecraft’s  perigee  into  the 
atmosphere  in  order  to  retard  its  velocity  and  shorten  its  semimajor  axis,  will  allow 
the  SMV  to  undertake  geotransfer  missions  and  make  it  back  with  enough  fuel  to 
reenter  and  land.  In  order  to  aerobrake,  the  SMV  would  need  to  perform  a  small 
burn  at  apogee,  AV  =  0.038  km/s,  in  order  to  dip  its  perigee  into  the  atmosphere. 
Once  the  SMV  finished  its  aerobraking  maneuver,  it  would,  if  needed,  perform  a  small 
burn  to  raise  its  perigee  above  any  noticeable  atmospheric  effects.  After  completing 
the  aerobrake  maneuver  and  circularization  burn,  it  would  be  in  a  near  circular,  low 
earth,  phasing  orbit.  The  SMV  could  proceed  directly  from  its  phasing  orbit  into  a 
reentry-to-land  flight  profile  or  could  orbit  the  earth  several  times  before  reentering 
in  order  to  minimize  the  required  maneuvering  in  the  continuum  flight  regime. 

Research  Objectives 

The  purpose  of  this  study  is  to  explore  the  possibility  that  the  SMV  can  fly 
out  to  a  target  in  geosynchronous  orbit,  perform  a  set  of  preprogrammed  tasks  on 
its  target,  and  return  to  a  circular,  low  earth  orbit  in  no  more  than  two  orbits  by 
aerobraking.  This  simulation  assumes  that  the  SMV  has  travelled  to  its  target,  suc¬ 
cessfully  completed  its  task,  and  now  needs  to  return  to  low  earth  orbit.  Aerobraking 

GEOSYNCHRONOUS 


Figure  9:  One  and  two-pass  aerobraking  orbits  [15,  55] 
will  allow  the  SMV  to  expend  a  majority  of  its  fuel  on  the  outbound  trip  instead 
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of  having  to  save  a  great  deal  of  it  for  the  return  trip  and  landing  procedure.  This 
greatly  enhances  its  operational  capabilities  and  makes  finding  a  viable  orbit  tra¬ 
jectory  of  the  utmost  importance.  Figure  9  depicts  two  potential  aerobraking  orbits. 
There  are  many  possible  aerobraking  schemes,  but  this  study  will  focus  on  two,  ballis¬ 
tic  entry  and  the  “Double-Dip”  reentry  method  described  in  Chapter  Eight  of  Wiesehs 
Spaceflight  Dynamics.  This  study  hopes  to  find  the  window  of  perigee  altitudes  that 
allows  the  SMV  to  maneuver  from  geotransfer  orbit  to  a  low  earth  circular  orbit  in 
no  more  than  two  passes  using  both  reentry  algorithms. 

Chapter  II  will  deal  extensively  with  the  problem  setup.  It  begins  by  explaining 
the  simulation’s  atmosphere  model,  its  planar  earth  assumption,  the  SMV  model,  the 
SMV’s  aerodynamic  properties,  and  the  mechanics  of  the  simulation’s  numerical  inte¬ 
grator.  The  development  of  the  equations  of  motion  for  both  reentry  schemes  follows 
the  discussion  of  the  assumptions.  Chapter  II  concludes  with  a  detailed  explanation 
of  the  “Double-Dip”  reentry. 

Chapter  III  delves  into  the  mechanics  of  the  simulation.  It  begins  by  explaining 
how  the  simulation’s  initial  and  final  conditions  were  determined.  It  will  also  detail 
the  algorithm  used  to  build  the  computer  simulation.  Chapter  III  ends  with  discussion 
of  each  subroutine. 

Chapter  IV  will  detail  the  simulation’s  results.  It  begins  by  discussing  the 
determination  of  the  bifurcation  point  using  both  ballistic  entry  and  “Double-Dip” 
reentry.  It  will  then  discuss  the  successful  perigee  altitude  choices,  disregarding  any 
constraints,  for  the  ballistic  entry.  Potential  “Double-Dip”  solutions  will  be  reviewed. 
Immediately  following  the  “Double-Dip”  solution  discussion,  the  required  technologi¬ 
cal  advancements  for  both  solutions  to  work  will  be  discussed.  Plots  of  the  experience 
accelerations  and  the  heating  rate  will  be  shown  to  display  any  SMV  shortcomings. 
Plots  of  the  SMV’s  flight  path,  velocity,  and  altitude  will  be  given  in  order  to  verify 
successful  mission  completion.  Chapter  IV  will  end  with  successful  perigee  window 
atmospheric  sensitivity  analysis. 
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Chapter  V  will  summarize  the  results  of  applying  a  ballistic  and  “Double-Dip” 
reentry  scheme  to  the  SMV’s  return  to  low  earth  orbit.  The  results  are  interpreted 
to  display  the  need  for  an  active  control  system  on  board  any  vehicle  performing  an 
aerobraking  maneuver  with  a  short  allowed  mission  completion  time.  Chapter  V  ends 
by  noting  areas  of  model  improvement  and  future  areas  of  study. 
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II.  Problem  Setup 


Chapter  Overview 

Chapter  II  starts  by  discussing  the  physical  assumptions  made  to  simplify  the 
problem  and  goes  on  to  discuss  the  SMV’s  aerodynamic  characteristics.  After  estab¬ 
lishing  the  framework  for  a  simulation,  the  simulation’s  numerical  integrator  is  iden¬ 
tified  and  discussed.  It  is  immediately  followed  by  the  development  of  the  equations 
of  motion.  Finally,  the  details  of  ballistic  and  “Double-Dip”  reentry  are  discussed. 

Physical  Assumptions 

It  is  important  to  establish  and  justify  the  simulation’s  assumptions  concern¬ 
ing  the  atmosphere,  earth’s  geometry,  aerodynamic  forces,  and  SMV  characteristics 
before  constructing  the  model  because  they  greatly  impact  the  performance  and  ac¬ 
curacy  of  the  model.  While  this  thesis  is  a  proof-of-concept  simulation  and  does  not 
need  to  include  every  minute  perturbation,  future  models  used  to  validate  operational 
capabilities  would  be  wise  to  include  third-body  gravitational  forces,  solar  wind  pres¬ 
sure,  the  nonspherical  earth  gravitational  accelerations,  and  any  vehicle  thrust.  These 
perturbations  are  extremely  small,  but  they  would  need  to  be  added  to  the  simulation 
in  an  application  that  requires  a  high  degree  of  accuracy. 

Atmospheric  Assumptions.  This  simulation  depends  on  the  interaction  be¬ 
tween  lift  and  drag  to  place  the  SMV  in  a  low  earth  circular  orbit.  Nothing  is  more 
important  in  determining  lift  and  drag  than  atmospheric  density,  and  nothing  could 
be  more  difficult  to  accurately  discern.  The  atmosphere  is  highly  volatile,  and  the 
density  at  any  location  is  affected  by  latitudinal  and  longitudinal  variations,  fluctu¬ 
ations  in  solar  activity,  magnetic  storms,  and  ocean  tides,  etc.  [22,  524-526].  The 
simulation  was  run  with  a  combined  Vallado  and  Regan  atmosphere.  To  highlight 
the  need  for  an  active  control  system  on  this  mission  because  of  atmospheric  fluctua¬ 
tions,  any  successful  orbits  were  integrated  again  assuming  a  five  percent  increase  in 
atmospheric  density.  The  sensitivity  analysis  will  be  discussed  in  Chapter  IV. 
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The  simulation  uses  a  combination  of  the  layered  exponential  atmosphere  de¬ 
scribed  by  Vallado  and  the  exponential  atmosphere  described  by  Regan  and  Anan- 
dakrishnan.  It  uses  Regan’s  and  Anandakrishnan’s  model  below  120  km  and  Vallado’s 
above  120  km.  Below  120  kilometers  [20,  38] 


P 


P  o  e 


h/H 


(2) 


where 

Po  =  Sea  level  density  (1.752  kg/m3) 
h  =  SMV  altitude  at  any  time  below  120  km  (km) 

H  =  Scale  Height  (6.7  x  103  m) 

At  altitudes  above  120  km,  Vallado’s  model  breaks  the  atmosphere  into  smaller 
layers,  each  with  a  unique  nominal  density  and  scale  height.  Above  120  km,  the 
atmospheric  density  is  found  by  using  the  following  equation  as  a  model  for  each 
altitude  layer  [22,  532-534]: 

/  heiip  Hq  .  .  . 

P  =  Po*  exp( - — - )  (3) 


where 

p  =  Atmospheric  density  at  altitude  (kg /km3) 

Po  =  Reference  density  at  h0  (kg /km3) 

ho  =  Reference  altitude  (km) 

heup  =  Height  above  the  ellipsoid  (km) 

H  =  Scale  height  at  altitude  layer  (km) 

Regan’s  atmosphere  model  closely  matches  the  1976  U.  S.  Standard  Atmosphere  up 
to  120  km  [20,  39].  Vallado’s  model  atmosphere  derives  its  values  from  25-500  km 
from  CIRA-72  and  uses  CIRA-72  with  exospheric  temperature,  T0 0  =  1000  K  for 
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500-1000  km  [22,  534],  Figure  10  shows  the  atmospheric  density  at  altitudes  from 
0-1000  km  using  the  combined  model. 


Figure  10:  Atmospheric  Density  Versus  Altitude  [22,  532-534] 

The  Matlab  code  used  for  the  simulation’s  combined  atmosphere  model  can  be 
found  on  the  compact  disc  that  accompanies  this  thesis.  It  is  also  assumed  that  the 
atmosphere’s  velocity  at  any  given  point  is  parallel  to  that  location’s  latitude  and 
does  not  act  in  the  radial  direction. 

Earth  Geometry.  The  earth  is  not  a  perfect  sphere,  but  the  perturbing 
accelerations  due  to  a  nonspherical  earth  are  negligible  compared  to  the  lift,  drag, 
and  gravity  forces  over  the  period  of  one  or  two  orbits.  Because  the  SMV’s  return 
from  geotransfer  orbit  will  only  span  one  or,  at  most,  two  orbits,  the  earth  is  modeled 
as  a  sphere  with  evenly  distributed  mass.  Future  operational  models  will  have  to 
account  for  these  accelerations. 

The  SMV’s  geotransfer  orbit  can  intersect  its  target  at  any  inclination,  and 
because  perturbing  forces  are  constrained  to  the  orbit  plane,  the  SMV’s  inclination, 
while  not  necessarily  zero,  can  be  ignored.  The  spherical  earth  assumption  allows  a 
further  simplification  to  a  two  dimensional,  circular  earth.  This  assumption  is  valid  as 
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long  as  the  lift  vector  acts  only  in  the  position-velocity  plane.  Conveniently  enough, 
this  simulation  constrains  the  lift  vector  to  the  position-velocity  plane. 


Aerodynamic  Forces.  The  question,  “When  should  atmospheric  perturba¬ 


tions  be  considered?”  is  difficult,  to  answer  because  the  addition  of  aerodynamic 
forces  increases  the  simulation  run  time  due  to  the  increasingly  complex  acceleration 
components  input  into  the  numerical  integrator.  However,  in  later  discussion  about 
the  “Double-Dip”  reentry  scenario,  atmospheric  entry  must  be  considered  well  above 
any  noticeable  accelerations.  In  order  to  maintain  validity,  lift  and  drag  are  considered 
when  either  of  their  magnitudes  are  greater  than  O.OOlg.  This  assumption  causes  the 
simulation  to  begin  considering  lift  and  drag  near  130  km.  The  off-on  “light  switch” 
method  reduces  the  numerical  integration  to  the  two-body  problem  when  lift  and  drag 
are  negligible. 

As  will  be  seen  below,  the  SMV’s  aerodynamic  characteristics  differ  in  each  of 
the  three  flowhelds:  continuum  flow,  transition  flow,  and  free  molecular  flow.  These 
flowhelds  are  defined  by  describing  the  relaxation  time  and  distance  for  their  con¬ 
stituent  molecules  [20,  314].  A  flowheld  with  a  high  number  of  molecules  has  a  high 
intermolecular  collision  frequency  and  short  relaxation  distance  and  time  [20,  313]. 
“If  this  relaxation  distance  is  small  compared  to  the  characteristic  dimension  of  the 
flowheld,  the  macroscopic  properties  can  be  considered  continuously”  [20,  313].  Free 
molecular  how  has  relaxation  distances  greater  than  the  characteristic  howheld  dimen¬ 
sion  [20,  314],  In  this  condition,  the  state  of  individual  molecules  must  be  considered. 
Thus,  Newtonian  impact  theory  is  the  governing  equation  for  this  howheld  [20,  314], 
A  nondimensional  parameter,  called  the  Knudsen  number,  was  created  to  help  differ¬ 
entiate  between  the  howhelds  [20,  314],  The  Knudsen  number  is  equal  to  [20,  314] 


characteristic  flow  field  dimension 


mean  free  path 


(4) 
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It  is  widely  held  that  Knudsen  numbers  >>  1  describe  free  molecular  flow,  and  Knud- 
sen  numbers  <<  1  pertain  to  continuum  [20,  316].  As  in  Captain  McNabb’s  study, 
“regions  where  Kn  >  10  will  be  considered  free  molecular  flow”  [13,  15],  and  “the 
region  where  Kn  <  .01,  the  flowheld  will  be  treated  as  continuum  flow”  [13,  15]. 
Transition  flow  resides  between  continuum  and  free  molecular  flow.  The  SMV  will 
traverse  all  three  regimes  as  it  dips  into  the  atmosphere.  Captain  McNabb  supplied 
equations  for  Cl  and  Cr>  in  molecular  free  flow  and  continuum.  In  order  not  to 
over  complicate  this  simple  model,  values  for  Cl  and  Cd  in  transition  flow  were  de¬ 
termined  by  linearly  interpolating  between  the  molecular  free  flow  and  continuum 
values.  Figure  11  presents  a  graph  of  altitude  versus  Knudsen  number. 


Figure  11:  Altitude  Versus  Knudsen  Number  [13,  16] 


SMV  Characteristics 

In  his  thesis,  Captain  Dennis  McNabb  dutifully  acknowledged  the  fact  that  the 
SMV  does  not  exist.  This  fact  sent  him  searching  for  potential  flight  vehicle  models. 
He  struck  engineering  gold  and  found  Boeing  SMV  concept  data  and  Boeing  X-37 
concept  drawings  [13,  11]. 

Captain  McNabb  used  concept  drawings  of  the  Boeing  X-37  to  deduce  its  body 
width,  nose  width,  and  chord  length  [13,  12],  He  then  used  these  figures  to  build  a  60- 
panel  model  of  the  SMV.  Appendix  B  shows  Captain  McNabb’s  60-panel  model.  This 
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Table  1:  Conceptual  SMV  Parameters  [13,  12] 


Parameter 

Value 

Weight 

44482.216  N 

Length 

8.8392  m 

Wingspan 

4.572  m 

Height 

2.8956  m 

Q  in 

X-37 

Length:  27.5  loot 

Wingspan:  15  feet 

Experiment  Bay:  7x4  loot 

1 

SI 

'JffrX  v 

Figure  12:  Boeing  X-37  [13,  12] 

model  was  used  to  determine  the  SMV’s  lift  and  drag  coefficients  at  varying  angles 
of  attack  and  flow  field  conditions.  He  performed  a  quadratic  regression  from  fifteen 
to  fifty  degrees  angle  of  attack  to  find  equations  for  the  lift  and  drag  coefficients. 
Appendix  C  contains  graphs  of  his  best  fit  lines.  His  equations  for  the  lift  and  drag 
coefficients  in  free  molecular  flow  and  continuum  are  given  below  [13,  20-22]: 


CLC  =  -3.3069  *  10 ~4(AOA)2  +  0.0552  *  (AO A)  -  0.5034  (5) 

CLF  =  -3.8241  *  10"5  *  (AO A)2  +  0.0033  *  (AO A)  -  0.0022  (6) 

CDC  =  9.6602  *  10"4  *  (AO A)2  -  0.0150  *  (AO A)  +  0.2402  (7) 

Cof  =  -1.3241  *  10"4  *  (AOA)2  +  0.0373  *  (AO A)  +  0.2265  (8) 


where 

Clc  —  Coefficient  of  lift  in  continuum 
Clf  —  Coefficient  of  lift  in  free  flow 
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Cdc  =  Coefficient  of  drag  in  continuum 
Cdf  —  Coefficient  of  drag  in  free  flow 

The  ballistic  simulation  used  a  Cdf  —  2.0  and  a  Cdc  —  3.75  in  continuum.  A  linear 
interpolation  supplied  Cd  values  for  the  transition  region.  Clc  and  Clf  =  0  for 
the  ballistic  reentry.  The  “Double-Dip”  reentry  simulation  was  run  twice  using  the 
SMV’s  best  Cl/Cd  angle  of  attack,  21°,  and  twice  at  the  maximum  Cd  in  Captain 
McNabb’s  angle  of  attack  window,  50°  [13,  22-23].  Appendix  C  displays  the  lift  and 
drag  coefficients  computed  from  McNabb’s  equations  and  their  ratios. 


Runge-Kutta  Method  of  Fourth  Order 


The  Runge-Kutta  Method  of  Fourth  Order  refers  to  numerical  integration  meth¬ 
ods  “originally  presented  by  Carl  Runge  (1856-1927)  in  1895,  and  Wilhelm  Kutta 
(1867-1944)  in  1901,  which  also  derive  from  a  Taylor  series”  [22,  500].  The  Runge- 
Kutta  Method  is  a  single-step  numerical  integrator;  this  means  it  “combines  the  state 
at  one  time  with  rates  at  several  other  times,  based  on  the  single-state  value  at  time, 
f0”  [22,  499].  It  calculates  those  rates  at  four  different  locations  in  each  time  interval 
using  a  fourth-order  Taylor  series  expansion  [22,  500].  “The  rates  are  readily  obtained 
from  the  equations  of  motion  and  allow  us  to  determine  the  state  at  succeeding  times, 
t0  +  /i”  [22,  499].  Before  the  Runge-Kutta  equations  are  given,  the  state  vector  must 


be  defined. 


where 


(9) 


f  =  SMV  inertial  position  vector 
v  =  SMV  inertial  velocity  vector 
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Taking  the  state  vector’s  derivative  with  respect  to  time  yields 


X  = 


(10) 


where 

a  =  SMV  inertial  acceleration  vector 

Vallado  defines  the  function  f(t,y)=X  =  y(t)  [22,  500].  The  Runge-Kutta  integrator 
will  calculate  the  values  of  X  using  the  following  equations  [22,  500]: 


2/i  —  /(h);  2/  o) 

h  h 

2/2  =  /Ro+2^o+2^i) 

r,hh 

2/3  =  /(h>  +  2 1 2/o  +  2^2)  (11) 

2/4  =  f(to  +  h,y0  +  hy3) 
yit )  =  y(t0)  +  \  (2/1  +  2j/2  +  2t/3  +  2/4)  +  0(h5) 

where 

h  =  Time  step  (s) 

0(h5)  =  Error  term  on  the  order  of  h 5 

It  is  immediately  obvious  that  time  step  selection  is  vital  to  producing  legitimate 
results.  A  test  orbit  without  perturbations  was  run  from  apogee  to  apogee  through 
the  Runge-Kutta  code  to  determine  the  most  appropriate  value  for  the  simulation’s 
time  step.  Table  2  shows  the  errors  associated  with  integrating  a  two-body  orbit  for 
one  orbit  with  five  different  time  steps. 

A  time  step  of  1  second  was  chosen  because  it  was  the  largest  time  step  to 
produce  less  than  0.5  nr  of  error.  Time  steps  larger  than  1  s  could  place  the  SMV  in 
the  wrong  atmosphere  layer  over  two  orbits  and  produce  erroneous  data.  The  error 
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Tabic  2:  Change  in  Position  and  Velocity  After  One  Orbit  at  Differing  Time  Steps 


Time  step 

0.25s 

0.5s 

Is 

2s 

5s 

A  X 

.00002  km 

.00010  km 

.00015  km 

.00052  km 

.00936  km 

A  Y 

.00006  km 

.00010  km 

.00018  km 

.00033  km 

.00145  km 

AA 

.00009  krn/s 

.00014  krn/s 

.00025  km/s 

.00048  km/s 

.00014  km/s 

associated  with  a  1  s  time  step  is  on  the  same  order  as  a  0.5  s  time  step.  Running  the 
simulation  with  a  time  step  smaller  than  1  s  would  lengthen  run  time  and,  considering 
the  simulation’s  assumptions,  not  necessarily  produce  better  results.  The  effects  of 
adding  perturbations  on  the  time  step  choice  must  also  be  considered.  Adding  pertur¬ 
bations  to  the  model  will  increase  the  model’s  error  because  the  perturbations  cause 
the  SMV’s  accelerations  to  change  rapidly.  However,  the  errors  in  the  SMV’s  position 
after  one  orbit  can  only  be  judged  relative  to  its  final  position  using  a  smaller  time 
step  because  the  SMV’s  true  position  and  velocity  after  one  orbit  can  never  be  known. 
Therefore,  position  and  velocity  errors  for  a  perturbed  orbit  were  not  considered  in 
the  time  step  selection.  Table  2  does  not  show  the  velocity  error  in  the  ey  direction 
because  it  is  smaller  than  0.5  mm  for  all  of  the  time  steps. 

Examination  of  Equation  11  reveals  that  the  SMV’s  acceleration  vector  or  equa¬ 
tions  of  motion  must  be  determined  before  the  simulation  can  begin.  This  will  be 
accomplished  in  the  next  section. 


Equations  of  Motion 

There  are  many  different  ways  to  formulate  the  equations  of  motion  for  the  SMV, 
and,  due  to  the  simulation’s  planar  earth  assumption  and  the  computing  power  of 
modern  personal  computers,  the  equations  are  formulated  using  Cartesian  coordinates 
in  the  geocentric-equatorial  coordinate  frame.  Hicks  describes  this  coordinate  frame: 

The  geocentric-equatorial  system  is  an  inertial  reference  frame  with  its 
origin  at  the  center  of  the  planet.  The  x-axis  points  in  the  direction  of 
the  vernal  equinox  and  the  z-axis  passes  out  the  north  pole.  The  y-axis 
completes  the  system  such  that  it  lies  in  the  planet’s  equatorial  plane  and 
&Z  &x  X  Cy.  [9] 
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This  coordinate  frame,  while  not  truly  stationary,  is  suitable  enough  to  be  considered 
inertial.  The  SMV  is  constrained  to  the  X-Y  plane  in  the  current  simulation  by 
cleverly  picking  its  starting  position  and  velocity  vectors,  but  the  simulation  has  the 
ability  to  upgrade  to  handle  position,  velocity,  and  acceleration  components  in  the  ez 
direction  by  retaining  zeros  in  that  position  that  can  later  be  changed. 

Newton’s  second  law  for  fixed  mass  sytems  will  be  used  to  find  the  SMV’s 
inertial  acceleration.  It  states,  “The  change  of  motion  is  proportional  to  the  motive 
force  impressed  and  is  made  in  the  direction  of  the  right  line  in  which  that  force  is 
impressed”  [2,  3].  Newton  was  simply  saying  that  the  force  on  an  object  is  equal  to 
the  time  rate  change  of  its  momentum,  or 


P  =  jtP  (12) 

Because  the  simulation  is  a  fixed  mass  system,  Newton’s  second  law  can  be  simplified 
into  the  more  recognizable 

F  =  ma  (13) 

where 

m  =  SMV  mass  (kg) 

Finding  the  SMV’s  acceleration  vector  now  becomes  a  matter  of  describing  the 
forces  acting  on  the  SMV  and  dividing  by  the  mass.  The  ballistic  reentry  simulation 
will  consider  gravity  and  drag,  while  the  “Double-Dip”  reentry  will  consider  gravity, 
drag,  and  lift. 

Gravity.  Newton  recognized  that  any  two  masses  act  on  each  other  with  a 
force  “proportional  to  the  product  of  their  masses  and  inversely  proportional  to  the 
square  of  the  distance  between  them”  [22,  135].  The  earth’s  mass  acts  on  the  SMV 
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along  the  SMV’s  inertial  position  vector  with  a  force  that  be  described  by  [22,  135] 


f9  =  - 


Gm^nrisat  f 


(14) 


where 

Fg  =  Force  of  gravity  (N) 

G  =  6.673  x  10”20  ±  0.001  x  10“20  ( km 3  •  kg~l  ■  s~ 2) 
mffi  =  5.9733320  x  1024  (kg) 

G  and  m®  are  often  combined  to  form  the  gravitational  parameter,  /x,  because  it  is 
easier  to  measure  the  gravitational  force  from  the  earth  felt  by  a  satellite  than  to  ac¬ 
curately  measure  the  mass  of  the  earth  and  the  gravitational  attraction  between  small 
objects  [22,  136].  The  simulation’s  computer  code  calculates  the  SMV’s  gravitational 
acceleration  with  the  following  equation: 


/i  r 


(15) 


The  simulation  was  run  taking  only  the  effects  of  gravity  on  the  SMV  into 
account  to  verify  its  construction.  The  SMV’s  orbit  should  be  one  complete  ellipse. 
Figure  13  validates  the  gravity  model. 


Lift.  Lift  plays  a  leading  role  in  the  ’’Double-Dip”  reentry  scheme.  The  SMV 
enters  the  earth’s  atmosphere  with  its  lift,  vector  pointed  away  from  the  earth,  and 
as  it  begins  to  pull  out  of  the  atmosphere,  its  lift  vector  is  rolled  to  point  towards 
the  earth.  This  roll  maneuver  is  usually  accomplished  at  a  pre-determined  flight  path 
angle.  There  are  multiple  ways  to  pick  the  flight  path  angle  at  which  the  SMV  will 
accomplish  its  “instantaneous”  roll;  these  will  be  discussed  later  in  the  chapter.  In 
keeping  with  earlier  assumptions,  the  SMV’s  lift  vector  is  constrained  to  the  orbit 
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Figure  13:  Two-body  Orbit 


plane.  The  resulting  simplified  planar  lift  is  described  by  the  equation  [23,  240]: 


(16) 


where 


Cl  =  Coefficient  of  lift 
A  =  Reference  area  (km2) 

Vrei  =  SMV  velocity  relative  to  the  atmosphere  (km/ s) 

The  force  due  to  lift  is  converted  into  an  acceleration  by  dividing  by  the  SMV’s  mass. 
The  acceleration  due  to  lift  can  be  expressed  by  the  equation: 


(17) 


2  m 


The  reference  area,  or  effective  lifting  surface,  was  determined  from  Captain  McNabb’s 
MATLAB  SMV  model.  Only  the  SMV’s  wings  were  considered  in  computing  the 
effective  lifting  surface,  but  future  studies  should  consider  including  the  wetted  area 
of  the  SMV’s  body  that  contributes  to  lift  production.  The  forward  lifting  surfaces 
were  broken  into  a  rectangle  and  triangle  to  more  accurately  determine  their  area. 
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The  two  lifting  surfaces  had  a  combined  area  of  approximately  3.875  m2.  This  lifting 
surface  area  was  used  for  the  21°  angle  of  attack  case  and  the  50°  angle  of  attack  case. 
Appendix  B  contains  pictures  of  Captain  McNabb’s  MATLAB  SMV  model. 

The  effects  of  a  roll  maneuver  with  lift  on  the  SMV’s  orbit  are  not  entirely 
intuitive.  Only  a  discussion  with  Dr.  Wiesel  validated  the  lift  perturbation  model.  The 
lift,  perturbation  rotates  the  orbit’s  apogee  clockwise  towards  its  original  perigee.  The 
lower  the  SMV  dips  into  the  atmosphere,  the  more  pronounced  this  effect  becomes. 
Figure  14  shows  the  effects  of  lift  on  the  SMV’s  orbit  when  its  perigee  dips  to  15  km, 
and  its  angle  of  attack  is  21°.  Figure  15  shows  the  effects  of  lift  on  the  SMV’s  orbit 
when  it  has  a  perigee  of  45  km  and  a  21°  angle  of  attack. 


x  io4  SMV  orbit 


Figure  14:  Lift  Perturbed  Orbit,  Perigee  =  15  km 


Drag. 

atmosphere. 


The  SMV  has  an  enormous  amount  of  kinetic  energy  as  it  reenters  the 
The  SMV  has  a  total  mechanical  energy  described  by  [23,  51]: 


£  = 


(18) 


and  a  kinetic  energy  per  unit  mass  T  =  v2  j  2  [23].  While  it  reenters  the  SMV  has 
a  specific  energy  of  roughly  -8.2  (km/ s)2  and  a  specific  kinetic  energy  of  50  million 
Joules  in  its  original  geotransfer  orbit.  In  a  400  kilometer  circular  orbit,  it  has  a 
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SMV  orbit 


x  104 


Figure  15:  Lift  Perturbed  Orbit,  Perigee  =  45  km 

mechanical  energy  of  approximately  -29  ( km/s )2  and  specific  kinetic  energy  of  30 
million  Joules.  This  means  the  SMV  must  find  a  means  to  discard  almost  21  (km/s)2 
in  one  orbit.  Drag  is  a  nonconservative  force  that  acts  in  the  direction  opposite  of  the 
SMV’s  velocity  that  will  bleed  away  the  SMV’s  energy  through  heat  transfer  to  the 
SMV  and  surrounding  air  molecules  [22,  521].  The  acceleration  due  to  planar  drag  is 
described  by  the  equation  [22,  522]: 

D  —  Cd  *  —  *  p*  Vre2  *  reZ.  (19) 

2  m  |  Vrel  I 

where 

Cd  =  Coefficient  of  drag 
S  =  Presented  area  (km2) 
nr  =  SMV  mass  (kg) 

The  amount  of  area  that  the  SMV  presents  to  the  oncoming  air  molecules  was  de¬ 
termined  using  a  flat  plate  model  for  each  of  the  lifting  surfaces  and  main  body.  At 
a  21°  angle  of  attack  the  SMV  presents  an  area  of  5.933  m2.  At  50°  angle  of  attack 
the  SMV  presents  an  area  of  12.682  m 2,  and  during  a  ballistic  trajectory,  the  SMV 
presents  16.555  m2 . 
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Drag  affects  the  SMV  by  shortening  its  semimajor  axis  and  reducing  its  eccen¬ 
tricity  [23,  85].  The  simulation  was  run  considering  only  drag  to  verify  the  correct 
construction  of  the  drag  model.  Figure  16  displays  the  SMV’s  orbit  around  the  earth 
and  verifies  the  drag  model. 


x  104  SMV  orbit 


km  x  104 

Figure  16:  Drag  Perturbed  Orbit 


Reentry  Schemes 

This  section  will  describe  the  two  reentry  techniques  used  in  this  study. 

Ballistic  Reentry.  When  the  SMV  flies  at  a  90°  angle  of  attack  it  has  effec¬ 
tively  zero  lift  and  plows  straight  through  the  atmosphere.  The  simulation  simply 
integrates  the  SMV’s  orbit  forward  from  the  target  satellite  and  determines  whether 
or  not  successful  final  conditions  are  met. 

“Double-Dip  Reentry”.  The  use  of  a  “Double-Dip”  reentry  dates  back  to  the 
American  Apollo  missions  and  Soviet  Zond  capsules.  The  Apollo  and  Zond  capsules’ 
offset  center  of  mass  created  a  lift,  vector  they  could  roll  to  maneuver  the  capsules 
in  the  atmosphere.  The  capsules  plowed  into  the  earth’s  atmosphere  with  their  lift, 
oriented  away  from  the  earth.  As  their  lift  began  to  raise  them  out  of  the  atmosphere, 
the  capsules  rolled  their  lift,  vectors  towards  the  earth.  Figure  17  displays  this  ma- 
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neuver  [23,  243].  In  Figure  17,  the  horizontal  axis  represents  downrange  distance,  and 
the  vertical  axis  represents  altitude.  Xr  and  Hr  denote  the  downrange  distance  and 
altitude  of  the  spacecraft  at  rollover. 


Down  Range  Distance 

Figure  17:  Potential  “Double-Dip”  Reentry  Flight  Profile  [23,  243] 

Had  the  capsules  not  altered  their  flight  profiles  after  their  initial  dip  into  the 
atmosphere  they  would  have  left  the  atmosphere  above  escape  speed  and  doomed  the 
astronauts  and  cosmonauts.  As  Figure  17  demonstrates,  the  solution  is  an  easy  one. 
It  was  necessary  to  roll  the  spacecrafts  to  point  their  lift  vectors  towards  earth  “to  help 
hold  the  vehicle  in  the  atmosphere  while  it  completes  its  deceleration”  [23,  243].  The 
capsules  performed  this  maneuver  during  the  upswing  portion  of  their  trajectory  to 
avoid  rocketing  the  spacecraft  into  the  ground  and  killing  any  of  their  highly  trained 
space  travellers.  The  roll  maneuver  is  an  easy  one  to  perform  and  has  been  thoroughly 
studied. 

Now  that  the  logic  behind  selecting  the  “Double-Dip”  reentry  scheme  has  been 
explained,  the  question,  “How  does  one  determine  the  flight  path  angle  at  which  to 
roll  the  SMV?”  must  be  answered.  Wiesel’s  explanation  of  how  to  choose  that  flight 
path  angle  will  be  reviewed.  Before  beginning,  it  is  important  to  define  the  flight  path 
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angle.  The  flight  path  angle  is  the  angle  between  the  SMV’s  velocity  vector  and  the 
local  horizon.  Velocity  vectors  residing  below  the  local  horizon  result  in  a  negative 
flight  path  angle.  It  is  necessary  to  switch  to  a  flat  earth  coordinate  system  at  this 
point  in  order  to  determine  the  flight  path  angle  at  which  to  roll  the  SMV. 

The  equations  of  motion  for  the  SMV  must  be  reconstructed  with  the  flat  earth 
coordinate  system.  The  reconstruction  will  begin  by  finding  the  force  of  gravity  on 
the  SMV.  “The  apparent  gravitational  force  on  the  vehicle  is  given  approximately 
by:”  [23,  239] 

F  =  _  mv 2  cos2 7 

^  ry*^2j  ry* 

As  one  can  see,  “the  apparent  force  of  gravity  would  be  negligible  compared  to  the 
aerodynamic  forces  on  the  vehicle  during  reentry”  [23,  239].  For  that  reason,  gravity 
will  be  ignored  for  this  formulation.  Now,  the  SMV’s  equations  of  motion  are:  [23,  239] 

dX 


dt 

dH 

dt 

dv 


m 


m  v 


dt 
d  7 

dt 


v  cosy 

(21) 

v  sin  7 

(22) 

—ad 

(23) 

L 

(24) 

where 


X  =  Downrange  Position,  (km) 

H  =  Altitude  (km) 

7  =  Flight  path  angle  (rad) 

Multiplying  the  above  equations  of  motion  by  [23,  240] 

dt  2m 


d'y  C^Apv 


(25) 


will  allow  them  to  be  integrated  over  the  new  independent  variable,  7.  As  long  as  the 
SMV’s  flight  path  angle  is  a  “monotonically  increasing  function  of  time”  [23,  240],  and 
it  is  assumed  to  be,  7  is  an  acceptable  independent  variable  and  the  SMV’s  equations 
of  motion  while  flying  right  side  up  are  [23,  239]: 


dX 
d  7 
dH 
d  7 
dV 
d  7 


cos  7  eH/H° 

Al 

(26) 

—  sin  7  eH/H° 

A  L 

(27) 

-—V 

Cl 

(28) 

where 


—  cl  Ap0 


2m 


Equation  28  can  be  rearranged,  using  separation  of  variables,  into  the  form  [23, 


240] 


dV  CD  , 

—  =  --cLil 


(29) 


Integrating  Equation  29  leads  to  [23,  240] 


fv  dV  CD  P  , 

lv.  ~  ~  ~C~L  L  7 


Vrsu  =  ve  exp  (  -^(7  -  7e) 


(30) 


where 

Ve  =  SMV  atmospheric  entry  speed 
7e  =  Atmospheric  entry  flight  path  angle 
Vrsu  =  SMV  right  side  up  velocity 

Equation  30  holds  true  until  the  SMV  performs  its  roll  maneuver  at  7r,  Hr ,  and 
Xr  [23,  243].  Realizing  that  the  only  change  to  the  equations  of  motion  is  the  direction 
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of  lift  allows  Equations  26  through  28  to  be  rewritten  as  [23,  243]: 


dX 
d  7 
dH 
d  7 
dV 
d  7 


1 

K~l 

1 


cos  7  e 


Ah 


szn  7  e 


cL 


H/H0 


(31) 

(32) 

(33) 


Substituting  a  velocity  at  rollover,  Vr,  as  the  initial  velocity  for  this  portion  of 
the  reentry  allows  Equation  33  to  be  integrated  like  Equation  28.  The  SMV’s  velocity 
while  upside  down  can  be  expressed  with  the  resulting  equation  [23,  98]: 

Vusd  =  K  exp  7  -  7r)^  (34) 

where 

Vusd  =  SMV  upside  down  velocity 
Solving  for  the  SMV’s  right  side  up  rollover  velocity 

vr  =  Ve  exp  ^-^(7r  -  7e)^  (35) 

presents  a  convenient  substitution 


Vusd  =  K  exp  (  -^E  (7r  _  %)j  x  exp  (  A 


(36) 


Solving  Equation  36  for  the  SMV’s  maximum  upside  down  altitude  results  in 


C 


Vusd{ 7  =  0)  =Ve  exp  (  ~^r( 2 7r  +  7e) 


(37) 
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A  targeted,  low  earth  orbit  with  an  altitude  of  300  km  would  provide  a  final  velocity 
of  about  7.7  km/s  per  second.  This  would  cause  Equation  37  to  resemble 


7.7  =  Ve  exp 


(  Ci(27r+7‘)) 


(38) 


Immediate  inspection  of  this  result  would  suggest  the  possibility  of  picking  a  final 
velocity  or  flight  path  angle  and  solving  for  the  remaining  variable  exists  because  7r 
is  a  function  of  7e.  Once  that  was  accomplished,  integrating  backwards  in  time  would 
provide  the  mandatory  initial  conditions;  however,  choosing  physically  feasible  final 
conditions  is  extremely  difficult,  and  for  that  reason,  this  method  was  not  tried. 


Taking  a  step  back  to  Equation  27  provides  a  simpler  approach.  Separating 
Equation  27  yields  [23,  241] 


(39) 


Integrating  Equation  39  gives  [23,  241] 


(40) 


where 


B  =  KLH0e~H^H°  -  cos'fe 

The  variable  B  is  introduced  at  this  point  to  simplify  the  notation.  As  Wiesel  points 
out,  “To  keep  the  argument  of  the  logarithm  greater  than  1  when  cos  7  =  1,  we 
must  require  that  1  +  B  <  KlH0”  [23,  241],  Solving  for  the  upside  down  altitude  by 
performing  a  similar  integration 


Husd  —  H0  In 


cosy  +  Busd 


(41) 
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To  ensure  that  Husd  is  greater  than  H0  the  natural  log  in  Equation  41  must  be 
positive.  This  means  [23,  244]: 

Busd  <  —  1  (42) 


Equation  40  and  Equation  41  are  equal  at  the  rollover  point.  Setting  them  equal 
yields  [23,  244] 

KlHq  =  -KlH0 

COS'Jr  +  BRsu  COS'Jr  +  BUSD 

Performing  some  simple  algebra  produces  an  equation  for  Busd  [23,  244]: 


Busd  =  —2  cos^r  —  BRSu 

Busd  =  cos^e  -  2cos^r  -  KLH0e~He/H°  (44) 


substituting  Equation  44  into  Equation  41  and  solving  for  the  flight  path  angle  at 
rollover  in  terms  of  the  desired  final  altitude,  H+  yields  [23,  244]: 


COS^fr  = 


1  +  cos'te  +  A'L.ffoO_"+/i'°  -  e-H‘/H°) 


(45) 


Equation  45  can  be  simplified  to 

1  . 

C0S7r  =  -  (1  +  COS7e) 


(46) 


because  H+  will  be  well  above  the  atmosphere  and  Hi  can  be  defined  above  the 
atmosphere. 

The  flight  path  angle  at  atmospheric  entry  is  below  five  degrees  for  every  case. 
The  small  flight  path  angle  at  entry  allows  the  two  cosine  terms  to  be  replaced  with 
the  first  two  terms  of  their  Taylor’s  series  [23,  252].  Equation  46  becomes 


1  - 


7r 


1 

2 


(47) 
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Simplifying  Equation  47  yields 


(48) 


and  because  ye  must  be  negative  for  reentry  to  occur,  and  is  positive  [23,  252]: 


7r  = 


1 

71 


7e 


(49) 


Equation  49  provides  a  convenient  way  to  find  the  approximate  flight  path  angle 
at  which  to  roll  the  SMV.  It  will  be  used  in  the  simulation  to  do  exactly  that.  Now 
that  an  expression  for  is  available,  it  can  be  substituted  back  into  the  expression 
for  upside  down  velocity  to  solve  for  the  flight  path  angle  at  atmospheric  entry.  For 
the  simulation  to  work,  the  final  upside  down  velocity,  Vusd  must  be  equal  to  the 
target  circular  velocity,  Vc .  Performing  the  substitution  yields  [23,  252]: 


Vc 


Ve  exp 


rearranging  Equation  50  gives  [23,  252]: 


CL  1 
CD  1  +  y/2 


In 


(50) 


(51) 


Equation  51  can  justify  this  study  if  it  proves  that  the  necessary  initial  condi¬ 
tions  are  physically  obtainable.  The  SMV’s  velocity  at  atmospheric  entry  is  roughly 
10-11.5  kilometers  per  second  while  the  target  final  velocity  is  below  7.75  kilometers 
per  second.  Inserting  those  values  into  Equation  51  gives  a  relatively  small  range  of 
negative  flight  path  angles  from  -13  to  -8.5.  This  answer  is  perfect;  the  flight  path 
angle  at  atmospheric  entry  must  be  negative  and  small,  and  it  is  both.  Only  numeri¬ 
cally  integrating  the  SMV’s  orbit  stands  in  the  way  of  finding  a  perigee  dip  altitude 
that  traps  the  SMV  in  a  circular,  low  earth  orbit. 
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Simulation  Parameters 


After  the  SMV  is  built,  it  will  have  parameters  that  cannot  be  violated  or  it  will 
be  destroyed.  This  study  will  propagate  each  orbit  without  regard  for  any  physical 
constraints  in  order  to  identify  areas  of  technology  that  need  improvement;  however, 
the  option  to  terminate  the  current  orbit  simulation  due  to  violations  of  the  theoretical 
SMV’s  constraints  exists  throughout  the  entire  simulation  and  can  easily  be  activated. 
The  simulation  tracks  the  SMV’s  acceleration,  heating  rate,  and  altitude  along  the 
entire  flight. 

Accelerations.  As  it  enters  the  atmosphere,  the  SMV  will  experience  an 
increased  acceleration  due  to  lift  and  drag,  but  it  does  little  good  to  know  just  the 
magnitudes  of  those  accelerations.  Most  aerodynamic  vehicles  express  their  maximum 
allowable  accelerations  in  terms  of  a  body  normal  acceleration  and  a  longitudinal 
acceleration.  The  Space  Shuttle  has  a  maximum  normal  load  of  2.5g  [13,  34].  This 
limit  will  not  be  imposed  upon  the  SMV  but  will  be  used  to  judge  the  SMV’s  body 
normal  acceleration  relative  to  the  Space  Shuttle. 

The  simulation’s  constant  angle  of  attack  assumption  allows  the  accelerations 
due  to  lift  and  drag  to  be  expressed  as  accelerations  in  the  body  normal  and  longi¬ 
tudinal  directions  by  performing  some  simple  trigonometry.  As  described  before,  lift, 
acts  perpendicular  to  the  vehicle’s  relative  velocity,  and  drag  acts  in  the  direction 
opposite  of  the  SMV’s  relative  velocity.  Figure  18  visually  displays  this  relationship. 


Figure  18:  Aerodynamic  Forces  on  Reentry  Vehicle  [13,  33] 
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In  Figure  18,  z  represents  the  body  normal  axis,  and  a  represents  the  longitu¬ 


dinal  axis  of  the  vehicle.  The  longitudinal  acceleration  becomes 


Longitudinal  Acceleration  =  Auft  sin(a)  —  Anrag  cos(a) 


(52) 


and  the  body  normal  acceleration  becomes 


Body  Normal  Acceleration  =  Auft  cos(ck)  +  Aorag  sin  (a) 


(53) 


where 

a  =  SMV  angle  of  attack 

The  body  normal  and  longitudinal  accelerations  felt  by  the  SMV  during  suc¬ 
cessful  orbits  will  be  displayed  in  Chapter  IV. 

Heating  Rate.  Heating  rates  during  atmospheric  reentry  are  one  of  the  largest 
design  constraints  for  any  reentry  vehicle.  The  SMV’s  heating  rate  could  be  quite  large 
because  it  enters  the  atmosphere  at  approximately  11  km/s.  Lu  provides  the  heating 
rate  constraint  equation  [12,  145] 


(54) 


where 


pref  =  1  kg/m3 
vref  =  1  m/s 
V  =  v/y/R0g0 

Qmax  —  Maximum  stagnation  point  heating  rate  ( W/m 2) 

Cq  =  Heat  flux  transmission  coefficient  {W/m2) 

To  maintain  continuity  with  Captain  McNabb’s  previous  research,  the  parameters 
for  maximum  stagnation  point  heating  rate  and  heat  flux  transmission  coefficient  are 
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taken  from  Lu  for  a  reusable  launch  vehicle  with  a  stagnation  point  heating  rate  based 
on  a  reference  sphere  with  a  radius  of  1  m  [12,  145]. 

qmax  =  544,300  W/m 2 

Cq  =  1.65  x  10_4  W/m2 

The  SMV’s  heating  rate  on  successful  orbits  will  be  discussed  in  Chapter  IV. 

Altitude.  Theoretically  there  is  no  limit  on  how  deep  into  the  atmosphere  the 
SMV  can  dip.  Once  it  is  built,  the  SMV  could  survive  any  number  of  reentry  profiles, 
but,  because  its  future  parameters  are  unknown,  two  altitude  floors  were  considered: 
45  km  and  70  km.  In  each  case,  the  simulation  was  ended  if  the  SMV  dipped  below 
the  altitude  floor.  Successful  dip  altitudes  will  be  discussed  in  Chapter  IV. 
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III.  Model  Construction 


Chapter  Overview 

The  purpose  of  this  chapter  is  to  explain  the  methodology  behind  the  simu¬ 
lation’s  construction.  It  will  detail  the  simulation’s  initial  and  final  conditions  and 
discuss  the  simulation’s  algorithm. 

Initial  and  Final  Conditions 

The  simulation  uses  apogee  of  the  SMV’s  geotransfer  orbit  as  its  initial  position 
because  the  SMV  will  return  immediately  after  completing  its  task.  The  SMV’s  initial 
velocity  is  found  by  targeting  a  perigee  altitude  and  substituting  it  into  Equations  57- 
59.  The  desired  final  conditions  were  formulated  keeping  the  goal  of  returning  the 
SMV  to  a  near  circular,  low  earth  orbit  in  mind.  Several  criteria  were  set  to  determine 
whether  or  not  the  current  orbit  was  successful.  To  be  considered  successful,  the 
SMV’s  altitude  had  to  be  between  184  km  and  600  km.  The  lower  altitude  was 
determined  from  investigating  the  lower  nominal  space  shuttle  mission  orbit  [19].  In 
addition  to  residing  within  an  altitude  window,  the  SMV’s  velocity  had  to  be  within 
500  m/s  of  circular  velocity  at  its  current  position,  and  in  order  to  ensure  the  SMV 
was  in  a  near  circular  orbit,  success  was  only  declared  if  its  flight  path  angle  was 
less  than  10°.  Orbits  from  the  ballistic  and  “Double-Dip”  reentries  that  meet  the 
above  criteria  will  be  discussed  in  Chapter  IV.  Any  target  perigee’s  orbit  that  did 
not  meet  the  success  criteria  after  being  integrated  through  two  atmosphere  passes 
was  considered  unsuccessful. 

Simulation  Algorithm 

The  work  required  to  build  a  robust  simulation  to  model  the  SMV’s  return  to 
low  Earth  orbit  is  not  finished,  and,  in  order  to  facilitate  future  work,  the  algorithm 
used  to  create  the  current  simulation’s  code  will  be  discussed  in  three  steps:  a  very 
broad  overview,  a  more  detailed  description,  and,  finally,  a  very  detailed  description 
of  each  function’s  code.  Figure  19  shows  the  simulation’s  algorithm. 
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The  simulation  is  comprised  of  one  main  function,  SMVdip,  and  four  subfunc¬ 
tions:  Cowelldip,  RK4dip,  Derivdip,  and  PDerivdip.  SMVdip  sets  up  the  simulation’s 
initial  conditions  and  calls  Cowelldip.  Cowelldip  contains  the  “Double-Dip”  reentry 
logic,  the  perigee  and  apogee  check  logic,  and  the  altitude  floor  check,  ft  also  calls 
RK4dip,  the  numerical  integrator.  RK4dip  calculates  a  new  state  vector  using  the 
returned  state  vector  time  derivatives  it  receives  from  Derivdip  and  PDerivdip.  De¬ 
rivdip  calculates  the  state  vector’s  time  derivative  without  perturbations.  PDerivdip 
calculates  the  state  vector’s  time  derivative  with  lift  and  drag  perturbations.  SMVdip, 
Cowelldip,  Derivdip,  PDerivdip,  and  RK4dip  can  be  found  on  the  compact  disc  that 
accompanies  this  thesis.  A  more  detailed  description  of  the  simulation’s  five  functions 
is  now  given. 

SMVdip  performs  several  important  operations  before  it  begins  integrating  the 
SMV’s  orbit,  ft  sets  the  simulation’s  start  date  and  time  to  00:00:00  UT  on  30  De¬ 
cember  2004.  The  time  and  date  were  arbitrarily  set,  but  future  simulations  that 
use  atmosphere  models  with  seasonal  variations  will  have  to  be  careful  in  setting  the 
simulation’s  time  and  date.  SMVdip  also  sets  the  integration  step  size  to  1  second. 
Drag  and  lift  perturbations  are  turned  off  or  on  in  the  simulation  depending  upon 
user  specified  choices  in  SMVdip.  The  circular  velocity  tolerance  and  starting  perigee 
altitude  are  also  assigned  values. 

At  this  point,  SMVdip  solves  for  the  orbit’s  initial  position,  velocity,  and  energy 
and  sends  them  to  Cowelldip.  Cowelldip  will  integrate  the  orbits  one  at  a  time  until 
it  finds  one  that  meets  the  desired  final  conditions.  Each  time  that  Cowelldip  returns 
an  unsuccessful  orbit,  SMVdip  increases  the  orbit’s  target  perigee  altitude  by  0.25 
kilometers,  recalculates  the  initial  position  and  velocity,  and  calls  on  Cowelldip  to 
integrate  the  new  orbit.  If  Cowelldip  returns  a  successful  orbit,  SMVdip  plots  the 
results  and  writes  the  solution  to  the  designated  output  hie. 

Cowelldip  begins  by  initializing  a  number  of  logic  variables.  After  those  variables 
are  set,  it  enters  a  loop  that  calls  RK4dip  until  the  desired  stop  time  is  reached.  Each 
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time  Cowelldip  calls  RK4,  RK4  returns  a  new  state  vector  one  time  step  after  the 
last  one.  Cowelldip  uses  the  new  state  vector  to  calculate  the  flight  path  angle,  and 
if  RK4dip  returns  a  variable  indicating  an  initial  encounter  with  the  atmosphere,  the 
flight  path  angle  is  stored  as  the  SMV’s  initial  flight  path  angle.  The  roll  flight  path 
angle  is  calculated  using  the  initial  flight  path  angle.  Once  the  SMV’s  flight  path 
angle  is  greater  than  the  roll  flight  path  angle,  Cowelldip  rolls  the  SMV  180  degrees. 
After  each  time  step,  Cowelldip  checks  for  perigee  and  apogee  passage,  Earth  impact, 
and  successful  final  conditions. 

RK4dip  needs  the  derivative  of  the  state  vector  with  respect  to  time  to  calculate 
the  new  state  vector.  RK4dip  calls  PDerivdip  and  Derivdip  to  find  the  derivative 
of  the  state  vector.  If  lift  and  drag  are  not  being  considered,  the  atmosphere  is 
“off” ,  and  Derivdip  creates  the  derivative  of  the  state  vector  using  only  gravitational 
acceleration.  Derivdip  still  calculates  lift  and  drag  but  only  uses  them  to  decide 
whether  or  not  to  turn  on  the  atmosphere.  PDerivdip  Ends  the  derivative  of  the  state 
vector  using  gravity,  lift,  and  drag.  PDerivdip  also  turns  “off”  the  atmosphere  if  the 
acceleration  due  to  lift,  or  drag  drops  below  the  preset  value.  Figure  19  gives  a  visual 
representation  of  the  simulation  algorithm. 

Detailed  Algorithm 

This  section  will  detail  each  step  in  the  five  MATLAB  functions  beginning  with 
SMVdip.  SMVdip  begins  by  converting  00:00:00  UT  on  30  December  2004  into  a 
Julian  date.  It  uses  a  function  called  JulianDay  written  by  Vallado.  The  JulianDay 
function  can  be  found  on  the  compact  disc  that  accompanies  this  thesis.  The  vari¬ 
ables  that  are  assigned  values  before  SMVdip  enters  its  loop  do  not  require  extensive 
calculations.  SMVdip  begins  its  loop  by  setting  up  the  orbit’s  initial  position  and  ve¬ 
locity.  Doing  this  requires  locating  its  target.  The  SMV’s  target  is  in  geosynchronous 
orbit  with  a  period  equal  to  twenty-three  hours,  fifty-six  minutes,  and  roughly  four 
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seconds.  The  target’s  orbital  radius  is  discerned  from  the  following  equation  [23,  56]: 


T  = 


(55) 


where 


T  =  Orbit  period  (s) 

H  =  398600.5  ( km3/s 2) 


To  simplify  the  simulation,  the  SMV’s  initial  position  vector  is  assumed  to 
coincide  with  its  target  at  its  apogee.  Constraining  the  SMV  to  its  orbital  plane 
creates  the  position  vector 


r  = 


42164.17 

0 

0 


km 


(56) 


SMVdip  then  reads  the  desired  perigee  altitude  to  dip  to,  and  the  simulation 
uses  that  information  to  place  the  SMV  in  an  equatorial  orbit  with  a  semimajor  axis 
described  by  [22,  104]: 

Ra  T  Ftp  /r>7\ 

a  =  - - -  (57) 

where 


Ra  =  Radius  of  apogee  (42164.17  km) 

Rp  =  Radius  of  perigee 

The  SMV’s  energy  can  then  be  found  by  substituting  the  semimajor  axis  into  the 
equation  [23,  56]: 


£  = 


k_ 

2  a 


(58) 


The  SMV’s  velocity  at  the  beginning  of  the  simulation  is  (Wiesel  1997:51) 


v  = 


(59) 
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and  resides  only  in  the  ey  direction  because  the  simulation  starts  at  apogee.  After 
calculating  the  SMV’s  velocity,  SMVdip  determines  the  orbit’s  end  time  with  the 
following  equation: 

End  time  =  Start  time  +  number  of  orbits  *  T/86400  (60) 


T  is  found  by  substituting  the  SMV’s  semimajor  axis  into  Equation  55,  and  the 
number  of  orbits  is  a  real  number  entered  by  the  user.  The  simulation  can  integrate 
forward  or  backward  in  time.  At  this  point,  SMVdip  calls  Cowelldip. 

Cowelldip  starts  by  initializing  logic  switches  and  then  immediately  calls  RK4dip. 
Cowelldip  tells  RK4dip  the  current  state  vector,  current  time,  which  perturbations  to 
consider,  and  the  direction  of  lift.  In  order  to  calculate  a  new  state  vector,  RK4dip 
calls  either  Derivdip  or  PDerivdip  four  times  at  three  different  time  steps:  once  at 
the  current  time,  twice  at  the  current  time  plus  half  the  step  size,  and  once  at  the 
current  time  plus  the  entire  step  size.  Both  subfunctions  calculate  the  accelerations 
due  to  gravity,  lift,  and  drag  at  the  SMV’s  current  position  and  velocity  to  formulate 
the  state  vector’s  time  derivative. 

PDerivdip  begins  by  calculating  the  gravitational  acceleration  with  the  equation 


X  = 


x 

y 


x  *  -If 


y*^f 
z  * 


(61) 


Once  PDerivdip  updates  the  state  vector’s  time  derivative,  it  must  find  lift  and  drag 
in  order  to  add  them  to  the  time  derivative.  The  SMV’s  velocity  relative  to  the 
atmosphere  must  be  found.  It  is  found  by  performing  a  vector  cross  product  with  the 
Earth’s  rotation  vector  and  the  SMV’s  position  vector.  Stated  in  equation  form 
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(62) 


where 


0 


u ;®  =  Earth’s  rotation  vector 


0 


r/s 


0.000072921151467 


The  SMV’s  velocity  relative  to  the  atmosphere  is  computed  with  the  equation: 


(63) 


The  next  step  in  computing  the  lift  and  drag  is  to  End  the  density  of  the  atmo¬ 
sphere.  PDerivdip  and  Derivdip  send  the  SMV’s  altitude  to  a  subroutine  called 
ATMOSPHERE2D  that  returns  the  atmospheric  density  at  that  altitude.  ATMO- 
SPHERE2D  can  be  found  on  the  compact  disc  that  accompanies  this  thesis.  Now  that 
the  SMV’s  velocity  relative  to  the  atmosphere  and  atmospheric  density  are  known, 
the  acceleration  due  to  drag  can  be  computed  using  Equation  19. 

Once  drag  is  computed,  PDerivdip  adds  it  to  the  state  vector’s  time  derivative. 

Now, 


x 


y 


(64) 


x  *  yr  +  Ad 
V  *  yr  +  Ad 


Derivdip  also  calculates  the  SMV’s  drag  but  does  not  add  it  to  the  state  vector’s 
time  derivative.  It  Ends  drag  only  to  determine  whether  or  not  to  turn  “on”  the 
atmosphere. 

PDerivdip  and  Derivdip  construct  the  lift  acceleration  vector  next.  Before  the 
SMV  rolls  over,  the  lift  vector’s  direction  is  determined  by  performing  a  vector  cross 
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product  with  the  SMV’s  relative  velocity  and  a  unit  vector  in  the  ez  direction.  The 
direction  of  the  SMV’s  lift  vector  is  then: 


L 


Vrel  X  k 


(65) 


where 


k  = 


0 

0 

1 


Multiplying  out  the  cross  product  reveals  that  the  lift  vector  points  away  from  the 
Earth.  L  is  turned  into  a  unit  vector  with  the  following  equation: 


L  = 


(66) 


Once  the  SMV  has  rolled  over,  Equation  65  becomes 


L 


Vrel  X  k 


(67) 


The  lift  direction  unit  vector  computation  does  not  change  from  above.  Once  the 
direction  of  lift  is  determined,  the  lift  magnitude  is  multiplied  by  the  lift  direction 
unit  vector  to  find  the  inertial  acceleration  experienced  by  the  SMV.  The  state  vector’s 
time  derivative  becomes 


x 


y 

z 

x  *  yr  +  Ad  +  Al 
V  *  yr  +  Ad  +  Al 
z  +  Ad  +  Al 


(68) 


43 


Like  before,  Derivdip  also  computes  the  SMV’s  lift  but  uses  it  only  to  determine 
whether  or  not  to  turn  “on”  the  atmosphere.  If  the  SMV’s  constraints  are  being 
considered,  PDerivdip  compares  the  magnitude  of  the  lift  and  drag  accelerations  to 
the  parameters  previously  discussed.  If  those  accelerations  surpass  the  limits,  flags 
in  PDerivdip  note  that  failure  and  force  the  simulation  to  start  over  with  the  next 
perigee  altitude.  Because  Derivdip  is  only  used  when  the  SMV  is  above  any  noticeable 
atmosphere,  it  does  not  contain  the  logic  flags  that  signal  constraint  violation. 

Derivdip’s  and  PDerivdip’s  final  task  is  to  calculate  the  heating  rate.  Equa¬ 
tion  54  is  simplified  to 

Heating  Rate  =  yfp  \  Vrei  |3  (69) 

If  SMV  constraints  are  being  considered  and  the  heating  rate  exceeds  the  maximum 
discussed  in  Chapter  II,  the  simulation  is  ended,  and  the  next  perigee  altitude  is 
attempted. 

PDerivdip  and  Derivdip  return  the  state  vector’s  time  derivative  at  three  differ¬ 
ent  times  to  RK4dip;  RK4dip  then  calculates  the  new  state  vector  using  the  equation: 

X2  =  X  +  (K1  +  2.0  *  (K2  +  K3)  +  Dt  *  A'4)/6.0  (70) 


where 

Dt  =  Integration  step  size 
K1  =  Dt*X(t0,X0) 

K2  =  Dt  *  X(t0  +  dt/ 2,  X0  +  Dt/ 2  *  Kl/Dt) 

K3  =  Dt  *  X(to  +  Dt/ 2,  X0  +  Dt/ 2  *  K2/Dt) 

K4  =  Dt  *  X (t0  +  Dt,  X0  +  Dt  *  K3) 

After  RK4dip  returns  the  new  state  vector,  Cowelldip  updates  the  time  and  finds 
the  SMV’s  flight  path  angle  by  calling  a  function  called  vecangle.  Vecangle  uses  the 
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equation  [11,  409]: 


cos  1  7  =  — — —  (71) 

|  b  ||  v  | 

where 

b  =  Vector  in  orbit  plane  perpendicular  to  SMV  radius  vector 

Cowelldip  corrects  the  sign  of  vecanglc’s  returned  flight  path  angle  by  multiplying 
it  with  a  logic  variable  that  equals  negative  one  while  the  SMV  travels  from  apogee 
to  perigee  and  positive  one  while  it  travels  from  perigee  to  apogee.  The  first  time 
Derivdip  detects  noticeable  accelerations  from  the  atmosphere,  Cowelldip  assigns  the 
current  flight  path  angle  value  as  the  atmospheric  entry  flight  path  angle.  The  roll 
over  flight  path  angle  is  determined  using  the  equation  [23,  252]: 


7 r 


1 

7! 


*  abs( 7e) 


(72) 


When  the  SMV’s  flight  path  angle  exceeds  the  roll  over  flight  path  angle,  Cowelldip 
rolls  the  SMV  by  telling  PDerivdip  to  use  Equation  67  in  place  of  Equation  65  when 
calculating  lift.  After  Cowelldip  checks  for  roll  over  it  determines  if  the  SMV’s  orbit 
meets  the  desired  final  conditions.  If  the  SMV  has  an  altitude  between  184  km  and 
600  km,  a  flight  path  angle  below  10  degrees,  and  a  velocity  less  than  500  m/s  smaller 
or  greater  than  circular  speed  at  its  current  altitude  that  orbit  was  successful.  Next 
Cowelldip  performs  an  apogee  and  perigee  check  and  an  Earth  impact  check.  If  the 
SMV  has  passed  perigee  twice  or  hit  the  Earth,  Cowelldip  stops  the  integration,  and 
SMVdip  picks  a  new  perigee  altitude  target.  Cowelldip  loops  until  success  is  acheived, 
the  SMV  hits  the  Earth,  or  the  SMV  passes  perigee  twice.  Cowelldip  reports  the  status 
of  the  current  orbit  to  SMVdip  once  it  is  finished  integrating.  If  the  pergiee  guess 
was  successful,  SMVdip  will  plot  the  SMV’s  orbit  around  the  Earth  and  update  the 
output  hie.  If  the  perigee  guess  was  unsuccessful,  SMVdip  will  increase  the  perigee 
altitude  by  250  meters  and  start  again. 
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Figure  19:  Simulation  Algorithm 
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IV.  Results 


Chapter  Overview 

This  chapter  will  present  the  results  obtained  from  numerically  integrating  the 
SMV’s  orbit  using  both  reentry  schemes.  Successful  one  and  two  pass  orbits  for  each 
reentry  scheme  will  be  examined.  The  effects  of  increased  atmospheric  density  on  the 
successful  perigee  windows  will  also  be  examined. 

Bifurcation 

Aerobraking  the  SMV  into  a  circular,  low  earth  orbit  in  one  or  two  passes  re¬ 
quires  dipping  deep  into  earth’s  atmosphere.  As  the  SMV  dips  its  perigee  deeper  into 
earth’s  atmosphere,  it  decelerates  faster  and  faster;  however,  this  increased  decelera¬ 
tion  does  not  come  without  risks.  If  the  SMV  dips  too  deeply  into  the  atmosphere, 
it  will  not  come  out.  The  bifurcation  point  is  the  lowest  altitude  that,  when  dipped 
to,  allows  the  SMV  to  exit  the  atmosphere.  If  the  SMV  were  to  dip  any  lower,  it 
would  be  forced  to  reenter  and  land  directly  from  geotransfer  orbit.  Because  this 
study  aims  to  return  the  SMV  into  a  low  earth,  phasing  orbit,  dipping  to  altitudes 
below  the  bifurcation  point  were  not  considered.  Identifying  the  bifurcation  point  for 
each  angle  of  attack  and  altitude  floor  identifies  the  beginning  of  potential  solution 
windows.  Table  3  shows  the  bifurcation  altitudes  that  were  found  using  empirical 
analysis  for  each  angle  of  attack  and  altitude  floor.  The  21°  and  50°  angle  of  attack 
bifurcation  points  were  found  using  the  “Double-Dip”  reentry  flight  profile.  The  90° 
angle  of  attack  results  represent  the  ballistic  reentry  flight  profile.  Two  altitude  floors, 
45  and  70  km,  were  used  in  the  simulation,  but  because  the  70  km  altitude  floor  did 
not  produce  any  solutions,  those  results  will  not  be  shown. 

Table  3:  Bifurcation  Altitudes  For  21°,  50°,  and  90°  Angle  of  Attack 


Altitude  Floor 

21° 

Or 

O 

o 

o 

O 

C7> 

45  km 

49.6  km 

48.9  km 

76  km 
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Successful  Simulations 

This  section  will  detail  the  successful  perigee  altitude  windows  for  each  reentry 
profile.  Because  no  successful  21°  angle  of  attack  perigee  altitudes  were  found,  21° 
angle  of  attack  orbits  will  no  longer  be  mentioned  in  this  chapter. 

Ballistic  Reentry.  Table  4  displays  the  successful  perigee  altitude  windows 
for  one  and  two  pass  ballistic  reentry.  In  Table  4,  the  target  altitudes  refer  to  the 
perigee  altitudes  used  to  formulate  the  SMV’s  unperturbed  velocity.  Once  the  SMV’s 
orbit  was  integrated  with  perturbations,  a  true  perigee  altitude  was  produced  for  each 
target  pergiee. 


Table  4:  Successful  Ballistic  Reentry  Altitude  Windows 


Number  of  Passes 

Target  Altitudes  (km) 

Actual  Altitudes  (km) 

1 

67.75  -  69.15 

66.801  -  68.449 

2 

72.55  73.75 

72.226  -  73.445 

The  lowest  and  highest  one  pass  solutions,  along  with  the  lowest  and  highest  two 
pass  solutions,  will  be  presented  because  they  most  likely  contain  the  lower  and  upper 
limits  for  the  SMV’s  heating  rate  and  accelerations.  Table  5  gives  the  SMV’s  altitude, 
instantaneous  velocity,  and  heating  rate  at  its  first  perigee  passage  for  each  altitude. 
Table  6  lists  the  SMV’s  experienced  accelerations  at  its  first  perigee. 

Table  5:  Instantaneous  Characteristics  at  Perigee  for  Ballistic  Reentry 


Altitude  (km) 

Passes 

Velocity  (km/s) 

Heating  Rate  (W/m2) 

66.801 

1 

9.240 

1005072.61 

68.449 

1 

9.506 

974556.91 

72.226 

2 

9.889 

832638.55 

73.445 

2 

9.969 

779789.93 

Table  6  gives  the  maximum  body  normal  acceleration  but  does  not  give  an  accurate 
representation  of  the  longitudinal  acceleration  experienced  by  the  SMV.  This  occurs 
because  the  acceleration  due  to  drag  acts  only  in  the  body  normal  direction  at  perigee. 
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Table  6:  Instantaneous  Accelerations  at  Perigee  for  Ballistic  Reentry 


Altitude  (km) 

Body  Normal  (g’s) 

Longitudinal  (g’s) 

66.801 

4.388 

0 

68.449 

3.651 

0 

72.226 

2.221 

0 

73.445 

1.863 

0 

Table  7  shows  the  approximate  maximum  value  for  longitudinal  acceleration  during 
the  maneuver  in  g’s. 

Table  7:  Approximate  Maximum  Longitudinal  Acceleration  for  Ballistic  Reentry 


Perigee  Altitude  (km) 

Passes 

Longitudinal  Acceleration  (g’s) 

66.801 

1 

2.75 

68.449 

1 

0.8 

72.226 

2 

1.4 

73.445 

2 

1.1 

It  is  also  important  to  analyze  the  successful  final  conditions  of  each  perigee’s 
orbit.  Table  8  details  the  SMV’s  final  altitude,  the  difference  between  its  velocity  and 
circular  velocity  at  its  current  altitude,  and  eccentricity. 


Table  8:  Final  Conditions  for  Ballistic  Reentry 


Altitude  (km) 

|  Vcirc  V actual  |  (kill/ s) 

Eccentricity 

One  pass  bottom 

184.11 

.0266 

.0342 

One  pass  top 

598.22 

.1971 

.2021 

Two  pass  bottom 

184.11 

.0328 

.0179 

Two  pass  top 

545.78 

.4999 

.1954 

These  results  show  that  as  the  SMV  dips  lower  into  the  atmosphere,  it  gets 
closer  to  leaving  the  atmosphere  in  a  circular  orbit  but  suffers  larger  accelerations 
and  heating  rates.  This  condition  causes  choosing  a  perigee  altitude  to  become  a 
trade  off  between  constructing  an  incredibly  strong  reentry  vehicle  that  needs  less 
fuel  to  complete  its  reentry  after  aerobraking  and  constructing  a  less  sturdy  vehicle 
that  must  retain  more  fuel  to  correct  its  reentry  velocity. 
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Figures  20-26  show  the  66.801  km  perigee  ballistic  reentry  orbit,  heating  rate, 
velocity,  altitude,  acceleration  in  the  body  normal  direction,  and  the  magnitude  of 
drag  through  the  entire  flight.  Figures  28-34  show  the  same  graphs  for  the  68.449  km 
perigee  ballistic  reentry.  Acceleration  due  to  lift,  in  the  ballistic  reentry  is  ignored 
because  it  equals  zero  for  the  entire  maneuver.  Figures  36  -43  show  the  same  graphs 
for  the  72.226  km  perigee,  two  pass  ballistic  reentry,  and  Figures  44-51  show  the  same 
graphs  for  the  73.445  km  perigee,  two  pass  ballistic  reentry.  All  of  the  following  figures 
can  be  found  in  Appendix  A. 

Investigation  of  Figures  20-26  show  that  the  entire  maneuver  is  accomplished 
in  less  than  six  hours.  Unfortunately,  Figure  24  reveals  that  the  SMV’s  heating 
rate  during  its  66.801  km  ballistic  pass  exceeds  the  qmax  established  in  Chapter  11 
by  a  factor  of  2.  Figure  32  shows  that  the  68.449  km  ballistic  pass  exceeds  the 
heating  rate  constraint  by  a  factor  of  roughly  1.8.  The  two  pass  solutions  fare  slightly 
better  by  only  violating  the  heating  constraint  by  a  factor  of  approximately  1.5.  If 
the  SMV’s  heating  rate  limit  resembles  Lu’s  constraint,  current  thermal  protection 
system  technology  will  not  suffice  in  returning  the  SMV  to  low  earth  orbit  in  one  or 
two  passes. 

“Double-Dip”  success.  The  “Double-Dip”  reentry  produced  successful  one 
and  two  pass  reentries.  Table  9  displays  the  successful  perigee  altitude  windows  for 
one  and  two  pass  “Double-Dip”  reentries.  In  Table  9,  the  target  altitudes  refer  to 
the  perigee  altitudes  used  to  formulate  the  SMV’s  unperturbed  velocity.  Once  the 
SMV’s  orbit  was  integrated  with  perturbations,  a  true  perigee  altitude  was  produced 
for  each  target  pergiee. 


Table  9:  Successful  “Double-Dip”  Reentry  Altitude  Windows 


Number  of  Passes 

Target  Altitudes  (km) 

Actual  Altitudes  (km) 

1 

73  -  77.85 

62.416  -  64.823 

2 

77.85  -  78.1 

64.823  -  64.962 
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The  lowest  one  pass  solution  and  the  highest  two  pass  solution,  along  with  the  highest 
one  pass  solution,  will  be  investigated  in  detail.  Table  10  gives  the  SMV’s  altitude, 
velocity,  and  heating  rate  at  its  first  perigee  passage  for  each  altitude. 


Table  10:  Perigee  Characteristics  for  “Double-Dip”  Reentry 


Altitude  (km) 

Passes 

Velocity  (km/s) 

Heating  Rate  (W/m2) 

62.416 

1 

9.240 

1413122.75 

64.8232 

1 

9.482 

1277526.36 

64.962 

2 

9.494 

1269114.22 

Table  11  lists  the  SMV’s  experienced  accelerations  at  its  first  perigee.  Displaying 
the  accelerations  due  to  lift  and  drag  will  help  determine  which  force  contributes  the 
greatest  to  the  SMV’s  acceleration. 


Table  11:  Perigee  Accelerations  for  “Double-Dip”  Reentry 


Altitude  (km) 

Body  Normal  (g’s) 

Longitudinal  (g’s) 

Lift  (g’s) 

Drag  (g’s) 

62.416 

3.028 

-1.549 

.764 

3.313 

64.823 

2.23 

1.142 

.5606 

2.446 

64.961 

2.19 

-1.120 

.550 

2.398 

Table  12  details  the  final  conditions  of  each  “Double-Dip”  orbit.  Figures  52-78  give 
the  SMV’s  orbit,  heating  rate,  velocity,  altitude,  acceleration  in  the  body  normal 
direction,  acceleration  in  the  longitudinal  body  direction,  the  magnitude  of  lift  ac¬ 
celeration,  and  the  magnitude  of  drag  acceleration  through  the  entire  flight  for  the 
62.416,  64.823,  and  64.962  km  perigee  “Double-Dip”  reentries.  Examination  of  those 
graphs  shows  that  for  the  lowest  possible  “Double-Dip”  reentry,  the  SMV  violates  the 
maximum  heating  rate  constraint  by  a  factor  of  roughly  2.75.  The  increased  heating 
rate  during  the  “Double-Dip”  reentry  is  caused  by  its  lower  dip  into  the  atmosphere. 
Stagnation  point  heating  analysis  would  reveal  an  even  greater  jump  in  the  heating 
rate  of  the  “Double-Dip”  reentry  compared  to  the  ballistic  reentry  because  of  its  re¬ 
duced  radius  of  curvature.  Figures  59,  68,  and  77  show  the  change  in  the  SMV’s  body 
normal  acceleration  as  its  perigee  is  raised  to  the  highest  two  pass  altitude. 
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Table  12:  Final  Conditions  for  “Double-Dip”  Reentry 


Altitude  (km) 

|  Vcirc  V actual  |  (km/s) 

Eccentricity 

One  pass  bottom 

184.042 

.0845 

.0520 

One  pass  top 

598.218 

.1930 

.1993 

Two  pass  top 

597.942 

.1934 

.1995 

The  results  from  this  study  can  be  compared  to  the  Powell  study  to  establish 
a  comparison  between  the  two  conceptual  reentry  vehicles.  Figure  79  reveals  that 
the  high  lift-to-drag  AOTV  must  dip  to  roughly  54  km,  8.4  km  deeper  into  the 
atmosphere  than  the  lowest  “Double-Dip”  reentry.  This  mission  requirement  most 
likely  stems  from  the  fact  that  the  AOTV  was  to  be  a  much  heavier  vehicle  than  the 
SMV.  It  had  to  dip  deeper  into  the  atmosphere  in  order  to  be  captured  in  a  low  earth 
orbit.  Dipping  deeper  into  the  atmosphere  should  cause  the  AOTV’s  heating  rate  to 
exceed  the  SMV’s  heating  rate.  Figure  81  shows,  once  the  units  are  converted,  that 
the  AOTV’s  heating  rate  does  indeed  surpass  the  SMV’s  envisioned  heating  rate  and 
exceeds  Lu’s  constraint  by  a  factor  of  11.  Figure  80  shows  that  two  vehicles  experience 
approximately  the  same  magnitude  of  accelerations. 

Atmospheric  Sensitivity 

The  ballistic  and  “Double-Dip”  simulations  were  run  a  second  time  with  a  five 
percent  increase  in  their  atmospheric  densities  to  demonstrate  the  perigee  windows’ 
susceptibility  to  atmospheric  fluctuations.  Table  13  shows  the  new  successful  perigee 
windows.  Increasing  the  atmospheric  density  by  five  percent  had  negative  effects  on 

Table  13:  Perigee  Windows  for  Increased  Atmospheric  Density 


Type  of  Reentry 

Number  of  Passes 

Altitudes  (km) 

Ballistic 

1 

67.046  -  68.513 

Ballistic 

2 

72.631  -  73.691 

“Double-Dip” 

1 

63.202  -  65.228 

“Double-Dip” 

2 

65.228  -  65.362 
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every  perigee  window  because  everyone  of  them  got  smaller.  Table  14  shows  these 
effects. 


Table  14:  Effects  of  Increased  Atmosphere  on  Perigee  Windows 


Type  of  Orbit 

Vertical  Change  (km) 

Width  Change  (km) 

1  Pass  Ballistic 

.245  Up 

Shrunk  by  .181 

2  Pass  Ballistic 

.405  Up 

Shrunk  by  .159 

1  Pass  “Double-Dip” 

.786  Up 

Shrunk  by  .381 

2  Pass  “Double-Dip” 

.405  Up 

Shrunk  by  .005 

The  crux  of  the  aerobraking  maneuver  is  the  SMV’s  burn  at  apogee  to  lower 
its  perigee  to  the  right  altitude.  As  the  perigee  windows  get  smaller,  performing  a 
burn  at  apogee  that  leads  to  a  successful  orbit  becomes  increasingly  difficult.  Table  15 
displays  the  differences  between  the  Ah’s  required  to  dip  to  the  bottom  of  the  perigee 
window  and  the  Aid’s  required  to  dip  to  the  top  of  the  perigee  window. 


Table  15:  AV  Requirements  with  Original  Atmospheric  Density 


Reentry  Type 

Perigee  Window  (km) 

A1 f bottom  ~  AVt0p  (km/ s) 

Ballistic 

66.801-68.449 

.000149 

Ballistic 

72.226-73.445 

.000122 

“Double-Dip” 

62.416-64.823 

.000516 

“Double-Dip” 

64.823-64.962 

.000027 

The  increased  atmosphere  shrinks  the  successful  perigee  window,  and,  as  a  re¬ 
sult,  increases  the  required  accuracy  of  the  burn  at  apogee.  Table  16  shows  the 
differences  between  the  Aid’s  required  to  dip  to  the  bottom  of  the  perigee  window 
and  the  Aid’s  required  to  dip  to  the  top  of  the  perigee  window  with  an  increased 
atmospheric  density. 

The  burn  at  apogee  to  enter  a  successful  orbit  with  a  known  atmospheric  density 
requires  a  great  deal  of  precision.  The  largest  Aid  window  spans  only  .5  m/s,  and 
the  smallest  Aid  window  ranges  only  .027  m/s.  Any  error  in  the  apogee  burn  greater 
than  the  AV  window  will  send  the  SMV  into  an  unsuccessful  orbit.  The  required 
Aid  accuracy  for  a  burn  with  uncertain  atmospheric  density  becomes  even  greater. 
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Table  16:  AV  Requirements  with  Increased  Atmospheric  Density 


Reentry  Type 

Perigee  Window  (km) 

bottom  A Vtop  (km/ s) 

Ballistic 

67.046-68.513 

.000133 

Ballistic 

72.631-73.691 

.000106 

“Double-Dip” 

63.202-65.228 

.000425 

“Double-Dip” 

65.228-65.362 

.000026 

The  change  in  both  perigee  windows  and  the  required  burn  accuracy  implies  that  a 
small  error  or  drastic  change  in  atmospheric  density  during  aerobraking  would  lead 
to  mission  failure.  At  this  point,  the  need  for  an  active  control  system  aboard  the 
SMV  becomes  glaringly  obvious. 
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V.  Conclusions  and  Recommendations 


The  purpose  of  this  study  was  to  prove  that  a  theoretical  SMV  could  maneuver  to 
geotransfer  orbit,  perform  a  task,  and  quickly  return  to  low  earth  orbit  in  no  more 
than  two  atmosphere  passes  with  its  advertised  AV  of  3.2  krn/s  by  aerobraking.  This 
was  done  by  numerically  integrating  the  SMV’s  orbit  considering  the  effects  of  lift, 
and  drag  around  a  planar  earth  model.  Now  that  this  task  has  been  completed,  the 
following  conclusions  can  be  drawn  from  the  results  provided  in  Chapter  IV.  This 
chapter  will  end  with  recommendations  for  future  study. 

Conclusions 

As  proven  in  Chapter  IV,  the  SMV  can  maneuver  to  geotransfer  orbit,  perform 
a  task,  and  swiftly  return  to  a  near  circular,  low  earth  orbit  using  aerobraking.  This 
ability  justifies  the  required  research  and  development  costs  because  it  will  greatly 
increase  the  SMV’s  mission  capabilities.  However,  the  SMV’s  ability  to  swiftly  return 
to  earth  with  a  minimal  amount  of  AV  comes  at  a  high  price.  The  SMV  experi¬ 
ences  extremely  large  heating  rates  that  will  require  a  thermal  protection  system  that 
outperforms  today’s  state-of-the-art  thermal  protection  systems,  and  the  rocket  burn 
needed  to  enter  a  successful  aerobraking  orbit  requires  an  accuracy  unavailable  from 
today’s  rocket  engines.  It  is,  however,  possible  to  overcome  the  hurdle  posed  by  the 
required  burn  accuracy  by  incorporating  a  number  of  smaller  correctional  burns  as 
part  of  an  active  control  system  as  the  SMV  approaches  perigee. 

The  ballistic  and  “Double-Dip”  reentries  have  significantly  different  characteris¬ 
tics.  The  “Double-Dip”  reentry’s  perigee  window  is  significantly  larger  than  the  bal¬ 
listic  reentry’s  perigee  window,  and  that  allows  for  greater  uncertainty  in  the  rocket 
burn,  but  it’s  greatest  heating  rate  is  larger  by  approximately  400,000  W/m2.  While 
the  “Double-Dip”  reentry  may  have  the  largest  perigee  window,  the  ballistic  reentry 
produces  the  orbit  with  the  best,  near  circular,  final  conditions.  The  lack  of  an  undis¬ 
puted  champion  between  the  two  reentries  makes  finding  a  solution  that  combines  the 
best  features  of  each  an  interesting  and  worthwhile  problem. 
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Recommendations 


1.  This  study  depended  on  keeping  the  SMV’s  angle  of  attack  constant.  Future  study 
of  this  problem  should  include  varying  the  angle  of  attack  as  a  function  of  velocity  or 
experienced  drag. 

2.  The  model  can  be  expanded  to  handle  three  dimensional  cases.  One  of  the  needed 
improvements  include  freeing  the  lift,  vector  from  the  position-velocity  plane.  Freeing 
the  lift  vector  from  the  position- velocity  plane  will  require  developing  expressions  for 
the  SMV’s  roll,  pitch,  and  yaw.  Changing  the  SMV’s  orientation  as  it  travels  through 
the  atmosphere  will  also  change  its  lift  and  drag  reference  areas.  Expressions  for  those 
must  also  be  developed.  Once  expressions  were  developed  for  the  SMV’s  three  con¬ 
trollable  angles,  it  would  be  more  advantageous  to  handle  the  reentry  as  a  two-point 
boundary  problem.  The  simulation  would  then  return  the  ideal  flight  trajectory  with 
detailed  profiles  of  its  roll,  pitch,  yaw,  and  angle  of  attack  during  the  flight. 

3.  Future  users  of  this  simulation  might  wish  to  include  smaller  perturbations  to  more 
accurately  model  the  SMV’s  orbit.  The  higher  order  geopotential  terms,  third-body 
gravitational  effects  from  the  sun  and  moon,  and  solar  radiation  are  just  a  few  of  the 
potential  perturbations  that  could  be  included. 

4.  Any  future  research  using  this  simulation  should  consider  working  in  a  different 
programming  language.  MATLAB  is  a  user  friendly  programming  language,  but  any 
advantages  in  programming  ease  are  lost  because  of  slow  computation  speed. 

5.  The  vehicle  parameters  used  in  this  study  are  all  estimates  for  the  SMV.  When 
a  prototype  is  built,  more  accurate  information  will  be  available  to  researchers.  Im¬ 
provements  in  the  SMV’s  aerodynamic  properties  should  be  included  as  they  are  made 
available  to  the  research  community. 
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Appendix  A. 


SMV  Orbit 


Figure  20:  Orbit  for  Ballistic  Reentry,  66.801  km  Perigee 

SMV  Orbit 


Figure  21:  Near  Earth  Portion  of  Orbit  for  Ballistic  Reentry,  66.801  km  Perigee 
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SMV  Velocity  Versus  Elapsed  Flight  Time 


Figure  22:  SMV  Velocity  Versus  Elapsed  Time  for  Ballistic  Reentry,  66.801  km 
Perigee 


Altitude  Versus  Elapsed  Flight  Time 


Figure  23:  SMV  Altitude  Versus  Elapsed  Time  for  Ballistic  Reentry,  66.801  km 
Perigee 
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Figure  24: 


Figure  25: 
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Figure  26: 


Figure  27: 


Acceleration  in  the  Body  Normal  Direction  Versus  Elapsed  Flight  Time 
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SMV  Orbit 


x  104 


Figure  28:  Orbit  for  Ballistic  Reentry,  68.449  km  Perigee 


SMV  Orbit 


Figure  29:  Near  Earth  Portion  of  Orbit  for  Ballistic  Reentry,  68.449  km  Perigee 
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SMV  Velocity  Versus  Elapsed  Flight  Time 


Figure  30:  SMV  Velocity  Versus  Elapsed  Time  for  Ballistic  Reentry,  68.449  km 
Perigee 


SMV  Altitude  Versus  Elapsed  Time 


Figure  31:  SMV  Altitude  Versus  Elapsed  Time  for  Ballistic  Reentry,  68.449  km 
Perigee 
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SMV  Heating  Rate  Versus  Elapsed  Time 


Figure  32:  SMV  Heating  Rate  Versus  Elapsed  Time  for  Ballistic  Reentry,  68.449 
km  Perigee 


Acceleration  Due  to  Drag  Versus  Elapsed  Time 
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Figure  33:  Acceleration  Due  to  Drag  Versus  Elapsed  Time  for  Ballistic  Reentry, 
68.449  km  Perigee 
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Figure  34: 


Figure  35: 


Acceleration  in  Body  Normal  Direction  Versus  Elapsed  Time 
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SMV  Orbit 
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Figure  36:  Orbit  for  Ballistic  Reentry,  72.226  km  Perigee 


SMV  Orbit 


Figure  37:  Near  Earth  Portion  of  Orbit  for  Ballistic  Reentry,  72.226  km  Perigee 
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SMV  Velocity  Versus  Elapsed  Flight  Time 


Figure  38:  SMV  Velocity  Versus  Elapsed  Time  for  Ballistic  Reentry,  72.226  km 
Perigee 


Altitude  Versus  Elapsed  Flight  Time 


Figure  39:  SMV  Altitude  Versus  Elapsed  Time  for  Ballistic  Reentry,  72.226  km 
Perigee 
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Figure  40: 


Figure  41: 


SMV  Heating  Rate  Versus  Elapsed  Time  for  Ballistic  Reentry,  72.226 
km  Perigee 
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Figure  42: 


Figure  43: 


Acceleration  in  the  Body  Normal  Direction  Versus  Elapsed  Flight  Time 
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SMV  Orbit 
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Figure  44:  Orbit  for  Ballistic  Reentry,  73.445  km  Perigee 


SMV  Orbit 


Figure  45:  Near  Earth  Portion  of  Orbit  for  Ballistic  Reentry,  73.445  km  Perigee 
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SMV  Velocity  Versus  Elapsed  Flight  Time 


Figure  46:  SMV  Velocity  Versus  Elapsed  Time  for  Ballistic  Reentry,  73.445  km 
Perigee 


Altitude  Versus  Elapsed  Flight  Time 


Figure  47:  SMV  Altitude  Versus  Elapsed  Time  for  Ballistic  Reentry,  73.445  km 
Perigee 
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Figure  48: 


Figure  49: 
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Figure  50: 


Figure  51: 


Acceleration  in  the  Body  Normal  Direction  Versus  Elapsed  Flight  Time 
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Figure  52:  Orbit  for  50°  AOA,  62.416  km  Perigee 


SMV  Orbit 


Figure  53:  Near  Earth  Portion  of  Orbit  for  50°  AOA,  62.416  km  Perigee 
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SMV  Velocity  Versus  Elapsed  Time 


Figure  54:  SMV  Velocity  Versus  Elapsed  Time  for  50°  AOA,  62.416  km  Perigee 


SMV  Altitude  Versus  Elapsed  Time 


Figure  55:  SMV  Altitude  Versus  Elapsed  Time  for  50°  AOA,  62.416  km  Perigee 
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Figure  56: 


Figure  57: 
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Figure  58: 


Figure  59: 
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Acceleration  in  Longitudinal  Body  Direction  Versus  Elapsed  Time 


Figure  60:  Acceleration  in  Longitudinal  Body  Direction  Versus  Elapsed  Time  for 
50°  AOA,  62.416  km  Perigee 


SMV  Orbit 


Figure  61:  Orbit  for  50°  AOA,  64.823  km  Perigee 
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SMV  Orbit 


Figure  62:  Near  Earth  Portion  of  Orbit  for  50°  AOA,  64.823  km  Perigee 


SMV  Velocity  Versus  Elapsed  Time 


Figure  63:  SMV  Velocity  Versus  Elapsed  Time  for  50°  AOA,  64.823  km  Perigee 
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SMV  Altitude  Versus  Elapsed  Time 


Figure  64:  SMV  Altitude  Versus  Elapsed  Time  for  50°  AOA,  64.823  km  Perigee 


Figure  65:  SMV  Heating  Rate  Versus  Elapsed  Time  for  50°  AOA,  64.823  km  Perigee 
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Figure  66: 


Figure  67: 
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Acceleration  in  Body  Normal  Direction  Versus  Elapsed  Time 


Figure  68:  Acceleration  in  Body  Normal  Direction  Versus  Elapsed  Time  for  50° 
AO  A,  64.823  km  Perigee 


Acceleration  in  Longitudinal  Body  Direction  Versus  Elapsed  Time 


Figure  69:  Acceleration  in  Longitudinal  Body  Direction  Versus  Elapsed  Time  for 
50°  AOA,  64.823  km  Perigee 
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Figure  70:  Orbit  for  50°  AOA,  64.961  km  Perigee 
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Figure  71:  Near  Earth  Portion  of  Orbit  for  50°  AOA,  64.961  km  Perigee 
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SMV  Velocity  Versus  Elapsed  Time 


Figure  72:  SMV  Velocity  Versus  Elapsed  Time  for  50°  AOA,  64.961  km  Perigee 


SMV  Altitude  Versus  Elapsed  Time 


Figure  73:  SMV  Altitude  Versus  Elapsed  Time  for  50°  AOA,  64.961  km  Perigee 


83 


Figure  74: 


Figure  75: 


SMV  Heating  Rate  Versus  Elapsed  Time 
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Figure  76: 


Figure  77: 


0.7 


Magnitude  of  Acceleration  Due  to  Lift  Versus  Elapsed  Time 


50  Deg  AOA 

64.961  km  Perigee 

- 

; 

- 

: 

- 

- 

,  1 

- 

0.5 


1  1.5 

Elapsed  Time  (s) 


2.5 


Acceleration  Due  to  Lift  Versus  Elapsed  Time  for  50°  AOA,  64.961  km 
Perigee 


Acceleration  in  Body  Normal  Direction  Versus  Elapsed  Time 


0.5 


50  Deg  AOA 
64.961  km  Perigee 


0.5 


1.5 

Elapsed  Time  (s) 


2.5 


Acceleration  in  Body  Normal  Direction  Versus  Elapsed  Time  for  50° 
AOA,  64.961  km  Perigee 


85 


Figure  78: 


Acceleration  in  Longitudinal  Body  Direction  Versus  Elapsed  Time 


± 


50  Deg  AOA 
64.961  km  Perigee 


1  1.5 

Elapsed  Time  (s) 


Acceleration  in  Longitudinal  Body  Direction  Versus  Elapsed  Time  for 
50°  AOA,  64.916  km  Perigee 


ALTITUDE  HISTORIES  FOR  MAXIMUM  RETURN  WEIGHT  AOTV’S 


Figure  79:  Altitude  History  for  Maximum  Return  Weight  AOTV  [15,  47] 


ACCELERATION  HISTORIES  FOR  MAXIMUM 


RETURN  WEIGHT  AOTV’S 


Figure  80:  Acceleration  Histories  for  Maximum  Return  Weight  AOTVs  [15,  48] 


REFERENCE  HEATING  RATE  HISTORIES  FOR 
MAXIMUM  RETURN  WEIGHT  AOTV’S 


Figure  81: 


Reference  Heating  Rate  Histories  for  Maximum  Return  Weight  AOTVs 
[15,  49] 
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Appendix  B. 


Measurements  are  scaled  to  the  vehicle  length:  L  =  29  ft  (8.84m) 


a.  0.1 74L 


b.  0.279L 


c.  0.220L 


d.  0.152L 


e.  0.175L 


Figure  50:  MATLAB  SMV  Model  Top  View  1 


a.  Nose  Length 

b.  Length  from  Nose  to  forward  lifting  surface 

c.  Forward  lifting  surface  inboard  length 

d.  Length  from  rear  of  forward  lifting  surface  to  front  of  rear  stabilizer 

e.  Length  from  front  of  rear  stabilizer  to  rear  of  vehicle 


Figure  82: 


MATLAB  SMV  Model  Top  View  [13,  73] 


Measurements  are  scaled  to  the  vehicle  length:  L  =  29  ft  (8.84m) 


f.  0.058L  g.  0.093L  h.  0.058L 


i.  0.238L 


Figure  51:  MATLAB  SMV  Model  Top  View  2 


f.  Rear  stabilizer  outboard  length 

g.  Rear  stabilizer  inboard  length 

h.  Forward  lifting  surface  outboard  length 

i.  Delta  wing  inboard  length 


Figure  83: 


MATLAB  SMV  Model  Top  View  2  [13,  74] 


Measurements  are  scaled  to  the  vehicle  length:  L  =  29  ft  (8.84m) 
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Figure  52:  MATLAB  SMV  Model  Front  View 


j.  Rotation  angle  for  Rear  Stabilizer 

k.  Nose  diameter 

l.  Body  width 

m.  Wingspan 


Figure  84:  MATLAB  SMV  Model  Front  View  [13,  75] 
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Appendix  C. 


Figure  85:  Coefficient  of  Drag  in  Continuum  vs  Angle  of  Attack 


Figure  86:  Coefficient  of  Drag  in  Free  Flow  vs  Angle  of  Attack 
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Figure  87:  Coefficient  of  Lift  in  Continuum  vs  Angle  of  Attack 


AOA  in  degrees 


Figure  88:  Coefficient  of  Lift  in  Free  Flow  vs  Angle  of  Attack 
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Table  17:  Lift  and  Drag  Coefficients  in  Continuum 


AOA  in  Deg 

cD 

cL 

CL/  CD 

1 

0.2261 

-0.4485 

-1.9831 

2 

0.2140 

-0.3943 

-1.8420 

3 

0.2038 

-0.3407 

-1.6713 

4 

0.1956 

-0.2878 

-1.4714 

5 

0.1893 

-0.2356 

-1.2446 

6 

0.1849 

-0.1841 

-0.9952 

7 

0.1825 

-0.1332 

-0.7297 

8 

0.1820 

-0.0829 

-0.4557 

9 

0.1834 

-0.0333 

-0.1819 

10 

0.1868 

0.0155 

0.0831 

11 

0.1920 

0.0637 

0.3320 

12 

0.1993 

0.1113 

0.5588 

13 

0.2084 

0.1583 

0.7594 

14 

0.2195 

0.2045 

0.9318 

15 

0.2325 

0.2501 

1.0758 

16 

0.2475 

0.2951 

1.1924 

17 

0.2643 

0.3394 

1.2838 

18 

0.2831 

0.3830 

1.3526 

19 

0.3039 

0.4260 

1.4016 

20 

0.3266 

0.4683 

1.4339 

21 

0.3512 

0.5099 

1.4520 

22 

0.3777 

0.5509 

1.4584 

23 

0.4062 

0.5912 

1.4555 

24 

0.4366 

0.6309 

1.4449 

25 

0.4689 

0.6699 

1.4285 
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Table  18:  Lift  and  Drag  Coefficients  in  Continuum 


AOA  in  Deg 

cD 

cL 

cL/  CD 

26 

0.5032 

0.7082 

1.4074 

27 

0.5394 

0.7459 

1.3828 

28 

0.5775 

0.7829 

1.3555 

29 

0.6176 

0.8192 

1.3265 

30 

0.6596 

0.8549 

1.2961 

31 

0.7035 

0.8900 

1.2650 

32 

0.7494 

0.9243 

1.2334 

33 

0.7971 

0.9580 

1.2018 

34 

0.8469 

0.9911 

1.1702 

35 

0.8985 

1.0235 

1.1390 

36 

0.9521 

1.0552 

1.1082 

37 

1.0076 

1.0862 

1.0780 

38 

1.0651 

1.1166 

1.0483 

39 

1.1245 

1.1464 

1.0194 

40 

1.1858 

1.1754 

0.9912 

41 

1.2490 

1.2039 

0.9638 

42 

1.3142 

1.2316 

0.9371 

43 

1.3813 

1.2587 

0.9112 

44 

1.4504 

1.2851 

0.8860 

45 

1.5213 

1.3109 

0.8616 

46 

1.5942 

1.3360 

0.8380 

47 

1.6691 

1.3605 

0.8150 

48 

1.7459 

1.3842 

0.7928 

49 

1.8246 

1.4074 

0.7713 

50 

1.9052 

1.4298 

0.7504 
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Table  19:  Lift  and  Drag  Coefficients  in  Free  Molecular  Flow 


AOA  in  Deg 

cD 

cL 

CL/  CD 

1 

0.2636 

0.0010 

0.0040 

2 

0.3005 

0.0042 

0.0141 

3 

0.3372 

0.0073 

0.0218 

4 

0.3735 

0.0103 

0.0278 

5 

0.4096 

0.0133 

0.0325 

6 

0.4455 

0.0162 

0.0364 

7 

0.4811 

0.0190 

0.0395 

8 

0.5164 

0.0217 

0.0421 

9 

0.5514 

0.0244 

0.0442 

10 

0.5862 

0.0269 

0.0460 

11 

0.6207 

0.0294 

0.0474 

12 

0.6550 

0.0318 

0.0486 

13 

0.6890 

0.0342 

0.0496 

14 

0.7227 

0.0365 

0.0505 

15 

0.7562 

0.0386 

0.0511 

16 

0.7894 

0.0408 

0.0516 

17 

0.8223 

0.0428 

0.0521 

18 

0.8549 

0.0448 

0.0524 

19 

0.8873 

0.0466 

0.0526 

20 

0.9195 

0.0485 

0.0527 

21 

0.9514 

0.0502 

0.0528 

22 

0.9830 

0.0518 

0.0527 

23 

1.0143 

0.0534 

0.0527 

24 

1.0454 

0.0549 

0.0525 

25 

1.0762 

0.0563 

0.0524 
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Table  20:  Lift  and  Drag  Coefficients  in  Free  Molecular  Flow 


AOA  in  Deg 

cD 

cL 

cL/  CD 

26 

1.1067 

0.0577 

0.0521 

27 

1.1370 

0.0590 

0.0519 

28 

1.1670 

0.0602 

0.0515 

29 

1.1968 

0.0613 

0.0512 

30 

1.2263 

0.0623 

0.0508 

31 

1.2556 

0.0634 

0.0505 

32 

1.2845 

0.0642 

0.0500 

33 

1.3132 

0.0651 

0.0495 

34 

1.3416 

0.0658 

0.0490 

35 

1.3698 

0.0665 

0.0485 

36 

1.3977 

0.0670 

0.0480 

37 

1.4253 

0.0675 

0.0474 

38 

1.4527 

0.0680 

0.0468 

39 

1.4798 

0.0683 

0.0462 

40 

1.5066 

0.0686 

0.0455 

41 

1.5332 

0.0688 

0.0449 

42 

1.5595 

0.0689 

0.0442 

43 

1.5856 

0.0690 

0.0435 

44 

1.6114 

0.0690 

0.0428 

45 

1.6369 

0.0689 

0.0421 

46 

1.6621 

0.0687 

0.0413 

47 

1.6871 

0.0684 

0.0406 

48 

1.7118 

0.0681 

0.0398 

49 

1.7363 

0.0677 

0.0390 

50 

1.7605 

0.0672 

0.0382 
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